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On Anomalous Variations of Maximum Electron Density 
and its Height in the F2 Region of the Ionosphere 


By Teruo SATO 


Geophysical Institute, Kyoto University, Kyoto 


Abstract 


An explanation of anomalous variations of the maximum 
electron density (proportinal to (f,F2)*) and its height (pF 2) of 
the F2 region is attempted. For this purpose, the reexamination 
of the behaviours of f,F 2, hpF2 and h'F2 in three seasons during 
sunspot maximum and minimum, is made and some interesting 
results are obtained. They are the followings. (i) The geomagne- 
tic control of /,F2 is conspicuous in low and middle latitudes 
in there seasons. But above middle latitudes the geographic 
control is remarkable. This is applied above 30° in winter 
(greatest secz) and equinox. (ii) fy#2 above 20° in December dis- 
tributes approximately in such a way as f)F 2<cos? x during sunspot 
minimum, but its distribution during sunspot maximum deviates 
from this line. (iii) The greater the zenith angle, the greater the 
increase of electron density from s.s. min. to s.s. max., and thus 
the magnitude of the seasonal anomaly of f,F2 is amplified far 
more during sunspot maximum than the minimum (iv) The-value 
of h,»F2 is greatest at geomagnetic equator throughout the year. 
(v) hpF2 becomes larger during sunspot maximum than the mini- 
mum. (vi) 4pF2 in December does not increase at higher latitudes 
and rather decreases. 

From these facts it seems that the variations of f,F2 and 
h,»F2 are controlled not only by the geomagnetic cause, but also 
by the cause which relates to the geographic latitude or solar 
activity. We investingate the geomagnetic distortion of daily varia- 
tion of electron density due to the vertical electron drift, taking 
the diurnal and semi-durnal components into consideration, and 
show that the trough in geomagnetic equator can be explained by 
this drift. The velocities of two kinds of drift are deduced directly 
from the data. A reason for geographic control is described. It is 
shown that the seasonal anomaly of f,f2 cannot be accounted for 
only by the vertical drift. Therefore in order to explain the 
seasonal anomaly of /,F2 and other phenomena, (iii)-(iv), mentioned 
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above, we suggest that there is at least the seasonal variation of 
the height distribution of the particle responsible for the F2 region. 
1. Introduction 

It has been. well known that the seasonal and daily variations of /,F2 (critical 
frequency of the F2 region) and h,F2 (height of the maximum electron density) 
of the F2 region) deviate remarkably from the law suggested by Chapmam [1]. 
One way to explain these anomalies is the thermal expantion of the region ((2], 
[3], [4]. However, the variation of the temperature to explain these anomalies is 
very great, therefore is not acceptable physically [5]. The variation due to the 
diffusion of an electron was considered by Ferraro [6]. According to his calculation 
the variation of the electron density due to this effect is very small even in the 
F2 region. The number density of particle taken by him is of two kinds. For 
smaller number density of particle the variation will become greater than that of 
his computation. But there still remains a doubt as to the behaviour of hpF2. 

Another excellent theory is the vertical electron drift by the electro-dynamic 
force firstly suggested and calculated by Martyn [7] and followed by Weiss [8] and 
Maeda [9]. This theory is qualitatively correct. It appears, however, that no 
satisfactory results have been found. 

Our recent study of the latitudinal distribution of f,F2 and h»F2 shows, 
besides the geomagnetic control of f,F2 which supported the drift theory, there is 
a geographic control in middle and higher latitudes, and the anomalies have a 
relation with the solar activity. Therefore we consider that the various kinds of 
anomalous variations are ascribed to the sum of the geomagnetic and the geo- 
graphic causes. It seems that at certain place and time the former is predominant ; 
and at another, the latter. In this paper, therefore, we first investigate whether 
or not the anomalies under the geomagnetic control can be accounted for by the 
vertical electron drift, and next in order to explain the anomalies which can not 
be explained by the electron drift, a suggestion is made. 

2. Statistical Results of Distributions of f,F2, h’F2 and h,F2 

For the purpose of well understanding and of later discussions we show in 
Fig. 1 (a)-(f) the latitudinal distribution of the value of f,F2, hpF2 and h'F2 (the 
minimum virtual height of the F2 region) at noon during the sunsopt maximum 
(s.s. max.) and minimum (s.s. min.). The upper two part in each fignre are respec- 
tively the distribution of f,F2 in geomagnetic and geographic latitudes, while the 
lower distributions are respectively the distribution of h’F2 (cross) and hpF2, 
cemuare} which are derived by Shimazaki’s [10] method, using the value of M(3000)F2. 
in ihn ees latitude. The data subjected to statistics are the hourly median 
values in ee Dec., 1947, Jun. 1948 and Jun., Sep., Dec. 1953. Jun., Sep., Dec. are the 
representatives of three seasons, the summer (May., Jun., Jul., Aug.), the equinox 
(Mar., Apr., Sep., Oct.) and the winter (Jan., Feb., Nov., Dec.). 
| It has been known that the value of foF2 shows less scatter when plotted 
in the geomagnetic latitude than in geographic ((11], [12], [13], [14]). But our 
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Fig. 1 (a)-(c) Distributions of noon median values of f)F2, h'F2 and hpF2 against geomagnetic and 
geographic (only /,F 2) latitndes. Period corresponds to sunspot maximum. 
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Fig. 1 (d)-(£) Distributions of noon median values of f,F2, h'F2, and hpF'2 against geomagnetic and 
geographic (only f,F2) latitudes. Period corresponds to sunspot minimum. 


detailed statistics show that though the above is consistent with our distribtion of 
foF2 in low and middle latitudes in geographic, the value of f,#2 above middle 
latitudes show better distributions in geographic than in geomagnetic. Especially, 
this is clearly seen in Dec. 1953 in northern hemisphere. The data in Dec. and 
Sep. 1952 and 1951 show the same result (Fig. 1 (A), (B), only for Dec.) There is 
a tendency that the smaller the average zenith angle at noon, of the month the 
better the geographic distribution. 

We have a reason to believe, as shown below, that in winter (greatest zenith 
angle) the value of f,F2 is well distributed in geographic latitude, even if the 
electron drift is present. Hence we plotted the value of f,F2 at noon at various 
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Fig. 1 (A) (B) Distributions of noon median value fp¥'2 against geomagnetic 
(uppes) and geographic latitudes in Decs. 1951 and 1952. 


stations agaist cosy in logarithmic scale. Fig. 2 shows f,F2-cosz% curve in Decs. 
1947, 1949, 1951-53. In these figures, a correction of grazing incidence is made [15]. 
It is evident that during s.s. mim. f)F2 increases linearly up to about 20°, and at 
low latitudes less than 20°, f,F2 decrease rapidly. This is the well known trough 
near magnetic equator and we must note that the peak of f,#2 near 20° indicates 
excessive electron density. 


The decrease of the electron density at the magnetic equator from the nor- 


i 


mal is estimated, from Fig. 2, as 


fifty percent. It is important 
that in general f,F2 in lower 
latitude deviates from the line 


foF2<cos? zy. The greater the 
sunspot number, the greater the 


deviation. 


[ , as usual relation, based on later 
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Fig. 2 Distribution of noon median value of f)F2 in discyagion: itoda; COnMIEGESS 
various latitudes (northern hemisphere) against that if the F2 region is formed 
cosy. Full line represents f)F'2cccos?y. by a different radiation from 

that of the #1 region and decay of the electron is of an attachment type, above 

relation is possible, but if the F2 region is formed in such a way as Bradbury’s 

[16] hypothesis, the above relation does not hold. This is because the production 

of the electron near the maximum of the F2 region in various latitudes is not so 

much different one another for the latter theory. However, above inferences 
are, needless to say, not decisive. 
There are also such notable characteristics as h,yF2 being greatest at geo- 


magnetic (or magnetic equator), and in general h»F2 increases far more during s.s. 
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max. than s.s. min. through two hemispheres, while /’F2 in geomagnetic equator 
decreases appreciablly. The harmonic analysis of the daily variation of h»pF2 at 
several stations shows that the augumentation of /,F2 is attributed to an incease 
of the constant trem (mean level) 40 
and not the sineand cosine terms. 

Next Fig. 3 shows the ratio of 30 
the electron density at noon of 


Sunspot Max. (1947-48) 


minimum coszx month to that of 20 
maximum. Above 23.5° this corres- 
ponds to the ratio of the electron 10 


density in winter to summr (opposite 
in souther hemisphere). During s.s. at 
min. this ratio is greater than untiy 
in middle and high latitudes in jqiueeees Sezer. ser 
northern hemisphere and at the rest 


of the stations less than unity. 


60 0 4 3 © 0 0 0 2 3% Bee Ce /) cy) 
Especially, it is interesting that in : Geographic Latitude 
higher latitudes than 60 this ratio is Fig. 3. Ratio (7) of noon median value of elect- 


ron density of the F2 region in minimum 
cosy month to that in maximum cosy 
month in various latitudes. 


less than unity. During s.s. max., 
the ratio increases world widely, 
except near the equatorial zone, and the augumentation of the ratio is greater for 
higher latitdes. This is due to the fact the increase of electron density from s.s. 
min. to s.s. max. is greatest in winter (smallest cosz) and smallest in summer 
(including equatorial zone) 
3. Some Results of Harmonic Analysis of h’F2..~ 

Harmonic analysis of h'F2 in summer and winter was made by Sato and 
Namikawa [17], using the data in about thirty observatories over the world. The 
phase and the amplitude of diurnal and semi-diurnal components are plotted agaist 
geographic and geomagetic latitudes. The data used there are the hourly median 
values of h'F2 in Jun., 1946 and Dec. 1945. The daily variation of h'F2 is analogous 
to that of h»F2, so that we can utilize these results. The results of analysis are 
summarized in the following. 

(i) In summer the amplitude of the diurnal component is greater in north- 


ern hemisphere than in southern and vice versa in winter. 
(ii) The amplitude of semi-diurnal component has an anologous distribution 


through two hemispheres as that in (i). But the difference of amplitude between 
two hemispheres at any solstice is smaller than that in the diurnal. 

(iii) The amplitude of diurnal component is greater than that of semi-diurnal 
in summer (winter) and is comparable with or slightly greater than that in winter 
(summer) in northern (southern) hemisphere. The magnitude of two components 
are respectively 50-70km and 30-40km in summer and 20-30km and 10-20km in 


winter. 
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(iv) The time of positive maximum of semi-diurnal component is within a 
range of one or two hours centered and at noon at midnight two hemispheres in 
two solstices. 

(v) In summer the time of positive maximum of diurnal component is within 
a range of two hours centered at noon in northern hemisphere and is within two 
hours after midnight in souhern. In winter these circumstances are approximately 
reversed. Near geomagnetic (magnetic) equator, the phase differs by 4-5 hours. 

The amplitude above mentioned is doubt, because the influence such as the 
retardation of wave is included. ‘The amplitudes of two components of hpF2 are 
smaller than those of h’F2, especially for high latitude. Near magnetic pole the 
amplitude of diurnal component h,F2 is 10-20km. 

4. Velocity of Electron Drift in the F2 Region 

Before advancing to a calculation of the variation of the electron densitiy 
due to the electron drift, we must deduce the velocity itself. Mitra [18] has com- 
puted this velocity from the night time data of the F2 region, assuming that the 
electron moves semi-diurnally. Judging from their results the velocity of the semi- 
drurnal drift is fairly small and is below 20km/hr. In this section we will obtain 
the velocity of not only the semi-diurnal but also the diurnal, by a method which 
is different from that of him. 

It we use the attachment law for the electron decay, the variaton of the 
electron density is given by 

oN =q—AN-09N, (1) 
where q is the electron production per c.c. per sec., v, the drift velocity of the 
electron which is independent of height, z, the coordinate of the height in unit of 
scale height. 8, strictly speaking, is not the attachment coefficient, but hereafter we 
call this as the attachment coefficient for the sake of convenience. We assume 
that the velocity of the electron is given by 

v=0,Sin (wi+A,)+v2Sin (2 wf+ 2), w=15°/hr. (23 


where v, and v, represent respectively the constant amplitude of two kinds of velo- 
city and 4,, 4, are phase lags. Substituting (2) into (1) and integrating the two sides 
of the equation from /, to /, and from /, to ¢;, we obtain two expressions. If the 
latter expression is subtracted from the former, the term concerning the electron 
production is eliminated when |f,—/.|= ¢—J/;! and we take f: as noon. Then we 
obtain the following expresson. 


ts ON es ON A ty ts to 
a—{ a= —4 | Nat—{ Nat|— (" ON ae egey 
\, ot dty ot J ty ty Jey ‘ dz vt AW a |. 

where fs is regarded independent of the height. Now we assume that the distribu- 


tion of the electron density during the integration is Chapman’s form, i.e 


N= js 3(1—z2—e7*) 
Nve : (4) 
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then 


Thus the expression (3) is writtien approximately as follows; 


3 N, 2 t p 9 
Ni — Ni = 4 bose = wae} (es ae 2 fet Nes ‘ : (I-2"), 3 


4 to cos (wt,;-+4,)+ 5 Vac08 (20f, + As)—0, COS (wh. + Ay) — ; Cos (2uts-+2s)! 


j 
— “Ahi e ; (1-e"*)ea. 91 14008 (ofs++ Ai) + : v»(2 wty-+ Ay) —v, COS (Wts+ A;) 


_ : v2CcOs(2 wt; +A2) i]. APs 


where the value at 4,2 is the mean of those at ¢, and #, and N,,, Ni. and N,,; are 
considered approximately equal. 

In oreder to deduce the velocity we adopt as the reference level the mini- 
mum virtual height of the F2 region. Since the critical frequency of the F1 region 
(f.F1) is about equal to the frequency corresponding to the electron density at 
that height of the F2 region (the difference between them is about 0.5mc/s), it 
appeears that the electron densitiy at h’F2 is equal to that corresponding to the 
frequency (f,F1+0.5)mc/s. At night the electron production is zero. Hence it is 
only sufficient to integrate the equation (1) from 4, to ¢, of any time. In that case 
the electron density at h’F2 is that corresponding to the frequency /,,;,,/ (minimum 
frequency reflected from the F region). 

It is noted that in progress of the calculation the electron density converted 
to a certain height should be used throughout the time of integration. The phase 
lags 4, and A. are obtained directly from the analysis of the same data from which 
the velocities v, and v, are calculated. 

The velocities of the diurnal and semi-diurnal drift at Kokubunji, Washington 
and Watheroo are shown in Table 1. At Kokubunji, two velocities are deduced 
from day and night time data. The data used here are the average of hourly 
median values of each month from Jan. 1950 to Dec. 1952. From the results in 
Table 1, it is found that v, is about 4-6m/s which almost equals to the value 


Table 1 
Kokubunji Washington Watheroo 
day . night day day 
V4 Vo V1 V2 Vy V2 Vy Do 
m/s! m/s m/s m/s m/s m/s m/s) m/s 
Summer Sol. eee Oo FO NS PLS 1.4 (33 
(May-Aug.) . ery eR eas ucla 
| Fon oe Sep., Oct.)| 14°] 52 | 13 | 54 | 18 | 4.2 | 1.7 | 3.2 | 
Winter. Sol. : ae 4.0 4,2 
(Jan., Feb., Nov., Dec.) oe | ae ost BG ae 
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obtained by Mitra. In general v, is about half of v. in winter and equinox at 
middle latitudes, while in summer v, is comparable with or greater than above 
and reaches to 6m/s. This fact shows that in summer the diurnal drift mainly 
contributes to the removal of the electron in middle and high latitudes. This state 
is consistent with the results obtained by harmonic analysis of h’F2. Near the 
geomagnetic equator, the behaviour of h’F2 is slightly different from that of hpF2. 
Hence above method cannot be applied. But judging from the variation of hpF 2, 
the diurnal velocity is probably greater than that in middle latitude and is 5-8m/s. 
5. Distortion of Maximum Electron Density Caused by Electron Drift 

In the previous section it was found that the diurnal drift is predominant 
over the semi-diurnal in summer (winter) in northern (southern) hemisphere. The 
distortion of the electron density so far investigated was for the semi-diurnal drift. 
In the present section we cmpute effect of two kinds of drift. The equation of 


the variation of the electron density is 

ON ON 
The electron production qg is ‘the function of the time and the height, and it is 
difficult to separate the two terms. For the convenience of the calculation we 
assume that 


(l—2-—e7*) 


php hae SP =q(a+bsinot)et—7-* ? (7) 


where q is the electron production at z=0 and x¥=0. a@ and } are constants which 
depend on the latitude. The assumption (7) implies that the production of the 
electron varies with coszx while the level of the maxmum production does not vary 
during the daytime. This invariance of the production level does not agree with 
the actual state which is usually repersented by Chapman’s law. According to 
the result of calculation in the case of no drift, variation of the electron density 
in the presence of drift under the condition of (7) appears to behave approximately 
as in the case of such a production as denoted by Chapman’s law. But the 
height variation at sunrise and sunset may differ from those in Chapman’s produc- 
tions, especially for the electron drift which begins at the same time as the electron 
production takes place, though such a difference is not much an important pro- 
blem here. Hence our calculations show a satisfactory result. 
We must also know the magnitude of @ at any height. It is considered that 
f can be obtained using the expression (6) by the method of least square. But the 
magnitude of f# thus obtained is very small and sometimes has the minus sign. 
So that the value is not reliable. Therefore we computed this value from the 
expression (= —~— N +, which is applied to the night time ionospheric data. WN is 
the maximum electron density when the height of the maximum does not much 
vary. The value of @ thus obtained is 1.5x10-*/s at Kokubunji; 2.2 10-°/s, at 
Washington ; and 1.7 10-*/s, at Watheroo at the level of 300km. 
Now since the velocity of the drift and the value of 8 are known, we advance 
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to compute the distortion of the electron density for the following two cases. 


(1) 2z>0 ~B=1x10-5 
z<0. - B=1x10-%e"= 
(2) z>0  B=2.72x10-5 (8) 


z<0 @=2.72x 10-e"* 
(i) Semi-Diurnal Drift 


The velocity of the drift is given by 
v=v,sin (2wt+2), (9) 


where 4 is the phase lag and w is constant. We divide the region into thin layer 
with constant density. Under such a condition the solution of (1) is given by 


4 282) sinot— wcos wt — 
N= gos Gaya) tL O(a) + Rew tO l™09 
Vo 


where = 2-5, L1— cos(2wi+)], Haagen S? ; (11) 
and @Q is the distributon of the electron density at ¢=0 (06500™). A(z), q(z), Q(z) 


implies that they are the function of z, and £(), q(~), Q(p) show that z is replaced by 
p. The daily variation of the electron density for two distributions of 6 and for a=0, 


1.0 


Local Time Local Time 


Fig. 4 Daily variations of electron density and Fig. 5 Daily variations of electron density and 
height of the moving region with uni- height of the moving region with uni- 
form velocity. Full line represents the form velocity. Full line represents the 
variation of the static region (no drift) variation of the static region (no drift) 
and R represents ratio of electron density and R represents ratio of electron density 
at any time to the maxiumum of the at any time to the maxmum of the static 
static region. 8=1%10-5/s for z>0, Ga region. $=2.72x10-5/s for z2>0; B=2,72 


«10-5e-7/s for 2<0. <10-5e-*/s for 2<0, 
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b=1 are shown in Figs. 4 and 5, in which the figures attached to the curves represent 
the time when the electron begins to move upwards or downwards. For example 
0t represents the electron drift which begins to move upwards at 00500™. The 
drift velocity used here is wo =4 which corresponds to 7.2m/s for H (scale height) 
=50km. The full line represents the variation in the case of no drift, which is 
called the static region. R is the ratio of electron density at any time to the 
maximum density of the day of the static region. 

It is clear from these figures that the smaller #, the greater the effect. It 
appears that the variation of height of maximum density obeys well to the move- 
ment of the electron. If we adopt far greater value of # than that used here, for 
example #=1x10-1/s, the effect of the drift is remarkably decreased and both the 
electron density and the height become closely similar to those of the static region. 

It is natural but important that the ratio of electron density at sunrise to 
the maximum of the day becomes greater for greater value of #8 and that the 
maximum density itself smaller. In the present, the maximum density in the static 


Diurnal 


0 2 4 2 OR OREO ie 


Local Teme 


4 16 18 2 22 24 


Local Time 


Fig. 6 Daily variation of electron density and Fig. 7 Daily variation of electron density and 
height of the moving region with uni- height of the moving region with uni- 


form velocity. Full line represents the 
variation of the sttaic region (no drift) 
and R represents ratio of electron density 
at any time to the maximum of static 
region. ~=1X10-5/s for 2>0; 
6=1X10-5e-*/s for 2<0. 


form velocity. Full line represents the 
variation of the static region (no drift) 
and #& represents ratio of electron density 
at any time to the maximum of static 
region. }=2.7210-5/s for 2>0; 
§=2.72x10-5e-*/s for 2<0. 
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region is respectively 3.8x10'q) and 1.93x104q, for cases (1) and (2). If the maxi- 
mum density of the F2 region is 1.5x10*/c.c. then q=37cm‘/s and 75cm?/s. The 
latter is approximatly equal to the electron production of the F2 region generally 
accepted (for example, Bates [19)). 

(ii) Diurnal Drift 

The velocity of the drift is given by v=wsin(wf+A). Then the variation of 
the electron density is represented by the same expression as (10). The distribution 
of f is taken as same as (8). The results are shown in Figs.6 and 7. It is evident 
drift gives far more influence than the semi-diurnal, and the electron drift of 6t 
is most effective. The reduction of the electron density relative to the static region 
reaches about fifty percent. 

According to the results mentioned in § 3, the diurnal compnent in summer 
and winter appears to correspond with the diurnal drift which begins respectvely 
to move up and down at midnight in northern himisphere. The variations of the 
height for 07 and 0) shown in Figs. 6 and 7 agree well with the diurnal components 
of h'F2. Therefore the electron may begin to move up and down at midnight in 
actual state. It may be considered, however, that as the radiation from the sun 
produces electron at high level at sunrise, the height for 6 rises more rapidly 
after the sunrise than that shown in Fig 7. In such a circumstance, the diurnal 
component of the daily variation of h'F2 may show also an apparent rise at mid- 
night. 

6. Election Drift in Atmosphere of Non-Uniform Temperature 

It is intereting to investigate an effect of the temperature with height gradient 

on the electron density when drift is occuring. Let the temperature increases 


linearly with height gradient 7. 
Then T= Ty +1 (h—y) (12) 


where T, is the temperature at the reference height ,. The variation of the 


electron density is given by 


aN aN _N OT 
ara tay ae (13) 


We use the similar approximate method as in § 5. For the diurnal drift, to which 
the velocity is given by v=wsin(ot+ 4), the solution of (13) is obtained as follows: 


ba 1% 
Zz caf be 5B(2) 7 Vo - eT 
N= fas (50+ 1 )eoswt— 9 eos “ot(* oe a Ta ) LQ(pje oT?) 


(% 
8 ha 1 flbelp) 1 _\)| FO Te) Cost é 
gat [( aE lis 5 (eet or ta (2)! for 4=0, (14) 


where p=2 —*1— —coswt), 
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and 
__9@) f A)[ (82 peel tind (282) 70% ) 
N=— om; COB LE m5 os | sin wt io +1 9 sin wt us oT (2) | 
ee i ot 1? LCL ae Oh ey Bee — ¢ sin‘ot) 
~ oT (2) (sinwt —5 sin ot)( 5 +1) 9 sin’ wf 6 sin’ wt 
( 156(2) , rv \]).-8t- ara Ar q(p)\ .—3(2)t— ore wt | 
- ae i) ale ais) |e + (»)+ w e 
for 4=+5, (15) 
where p=z- “a sin of. 
When A=z and A= =F , the solution is given by the same expressions as (14) 


and (15) except replacing v) by —v). The distributions of temperature and of coef- 
ficient 6 with height are shown in Table 2. The velocity adopted is5m/s. In this 
calculation it is supposed that the production maximum of the electron lies at 270km 
in summer (xz is smallest) and at 370km in winter. The electron drifts downwards 
in winter and upwards in summer. It is not evident whether or not this assump- 
tion agrees with the real condition, but the effect of the temperature will be 
verified. Figs. 8 and 9 show the results of calculations. The static regions in two 


Height 


Local Time 


Local Ti me 


Fig. 8 Daily variation of electron density and Fig. 9 Daily variations of electron density and 
height of moving region with uniform height of moving region with alee 
velocity when ~ temperature increases velocity when temperature increases 
linearly with height. Full line represnts linearly with height. Full line represents 
static region, the static region. 
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figures represented by the full lie are res- Table 2 

pectively similar with those in Fig. 6 and 7, Height Temperature| Attachment Coeff. 
and the maximum density of the day is 3.9 “400km| —-2200°k 1X10-5/s 
x10*qg and 2.3x10'g. The daily variation 380 | 2000 1.151075 

of the electron density in Fig 8 for 0) shows 360 ‘igoo. | )2=(1.25x10-5 tS” 
the great increase compared with the 340 | 1600 | 1.35x10-7— 
corresponding curve in Fig. 6. On the other oo WT eEO 1.5 X10-5 
hand, the variation of 67, 07, shows the “300 | 1200 1.8 x10-5 
decrease relative to the static region than 280 1000 2.2 10-5 

the corresponding curves in Fig 7, especial- b60- "1 B00” gO 

ly, the decrease for 67 reaches 50%. The “240 600 3.9 10-5 
variation of height is analogus to those in 990 | 400 6.0 x10-8 
Ga bund 7. 200 400 1.0 x10-4 


7. Explanation of Daily and Seasonal 
Anomalies of Electron Density and Height 
General aspects of variations of f,F2 and h’F2 are shown in Fig. 10.. Since 
the variations of f,F2 and h,F2 (analogaus to h’F 2) in winter in northern hemisphere 
are analogous to those in summer in the southern and the reverse case also holds, 
we hereafter disccuss the variations in the northern hemisphere. 
(i) We first investigate whether the daily variations in two solstices can be 
explained by the vertical drift. We consider that the effect of the semi-diurnal 
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Fig. 10 Daily variations of fF 2 and h'F2 in various latitudes in Jun. and Dec. 1953. 
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drift is small. As shown in § 5 and § 6, including the effect of the temperature 
gradient with height, the decrease of the electron density relative to that of the 
static region (no drift) is about 20% for the diurnal drift of 07 which is considered 
to occur in summer and about 50% for the drift of 6? which is cosidered to occur 
near magnetic equator, while for the drift of 0) which corresponds to the case in 
winter, the daily variation is considered to be almost equal to the static case. In 
equinox the drift is mainly of semi-diurnal and drift of 0) adds in high latitudes. 
So that in equinox the daily variation is approximately the same as the static 
region. But in this season the anomaly in electron density exists. This reason is 
inferred in (v) of this section. In winter since the daily variation, save the case 
of low latitude, is regarded as that of the static region, the observed variation will 
be obtained, provided that the attachment coefficient in the region of electron 
production is 3x10-°/s-5x10-°/s based on the estimation shown in § 8. On the 
other hand, in low latitudes throughout the year and above the middle latitudes in 
summer, the data indicate that an upward drift influences on the electron density and 
height. In § 2 we showed that near the magnetic equator the electron density 
is reduced by about 50% from the static Chapman norm. This is understood easily 
if we consider that there is an electron drift of 67. If this drift is present through- 
out the year, the trough in magnetic equator is always set up. Hirono and Maeda 
[20] also indicated that the vertical electron drift of 7.357 leads to an observed 
daily variation at Huancayo, and the reduction of the electron density from the 
static Chapman norm is abont is 40%-60%. Their method of integration of the 
equation of continuity is numerical and the decay of the electron is of a recombi- 
nation type. Our calculation differs more or less from those by them, but the results 
are well consistent with theirs. Now the summer daily variation of electron density 
become smaller for higher latitude, and the ratio of density at sunrise to the maximum 
of the day is 1/1.5-1/3, while 1/10-1/15 in winter. This makes us to regard that the 
static region in summer lies at the level of smaller value of § than that in winter. 
But in all our calculations the sunshine begins at 06500™ and sunset at 18h00™. In 
reality the time of sunrise approaches more to midnight as the latitude becomes higher, 
the duration of the sunshine, therefore, longer and the variation of the intensity 
of the radiation smaller. Under these conditions the daily variation of the electron 
density becomes smaller. As an extreme case we regard the sun shines all over 


Yet ot 
implies a constant electron density. This characteristic makes B to be small 


apparently. If this characteristic and the influence of the vertical electron drift of 
Of are taken into consideration the summer daily variation of the density above 


the day with a constant intensity. Then the solution ele a sen is © which 
B 


middle latitude can be explained when the attachment law is applied as the elect- 
ron decay. If the decay of the electron is recombination type, we cannot help but 
consider that the production of the electron in summer is higher than in winter. 
The height variation is accounted for by vertical drift. 
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(ii) Let the static region lies at the height respectively 4x 10-°/s and 6x 10-*/s 
in winter (w) and summer (s) about 30° in_ latitude (reason described in § 8). Then 
the ratio (3) of the electron density at noon in two solstices is * ip Gs. 
where gq is the electron production at noon. This estimation is based on the ae 
variation mentioned in (i). If the electron is proportional to cosy, N,=N, is 
obtained. At 50°, ¥~ 7 rm and since q,~34q,,, N,,<N\, is obtained. Thus the 
seasonal anomaly of the electron density in such latitudes does not occur. Never- 
theless, in reality, anomaly takes place. If we take Bradbury’s hypothesis, the maxi- 
mum of the electron density in the F2 region lies higher by two or more in scale 
height than that of the Fl region. If we take two in scale height, N,~1.2N,, 
at 50°, If distance of two maximum is more than two in scale height, the pos- 
sibility N,~WN,, is present. So that, a part of the seasonal anomaly is in this 
case explained. In consquence, the anomalous seasonal variation cannot be ac- 
counted for only by the vertical drift. 

(iii) There are other facts which can not be explained by the vertical drift 
only. They are: 

(a) Increase of electron density from s.s. min to s.s. max. is greater in writer than 


in summer and the seasonal anomaly of electron density is amplified world- 
widely during s.s. max. far more than s.s. min. as shown in Fig. 3. This ampli- 
fication therefore becomes larger for higher latitude. 

(b) There is an increase of h,F2 during s.s. max. compared with that during s.s. 
min.. 

(c) The latitudinal distribution of f,F2 even in equinox is better in geographic 
than in geomagnetic above 30° and a slight seasonal anomaly still exists. 

(d) hpF2 in winter does not incsease in higher latitudes and rather decreases. 

(e) Distribution of f,F2 in low latitude against cosy deviates from’the line fpF2 
ecos?z. The greater the sunspot number, the greater this deviation. 

(f) There are many stations where h'F1 in summer is larger than that in winter 
or at least equal, as shown in Fig. 11. This is the reverse to Chapma’s law. On 
the other hand, daily varia- 
tions of /F2 and f,F1 obey 
well his law. 

It is necessary to add 
brief comments. First to (a) and 
(b), we first considered that both 
the amplification of anomaly of 
election density and the increase 
of h»F2 are due to the augu- PP ars EON ENE NETS eres vt 
mentation of the vertical drift. Fig. 11 Seasonal variation of noon median value of h’F1 
Since it is likely that the vari- from Jan. 1951 to Dec. 1953. 
ation of hp»F2 represents a direct result of the drift, the daily variations of hpF2 
at Kokubunji, Washington and several other stations during s.s. max. and s.s. min. 
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are put in a harmonic analysis. The result shows that only the constant term 
(mean level) incseases and no appseciable increase of diurnal and semi-diuanal com- 
ponents takes place. To(c), it is known from a harmonic analysis of the height that 
there is no appreciable effect of the diunal drift in equinox. Hence the slight 
anomaly is probably due to the cause associated with the geographic latitude. 
Finally to (f), we consider that the behaviour of /#’/F1 in summer is due to an ele- 
vation of the production, because in such a height § is very large and the electron 
drift is remarkably reduced. 

At any rate in (ii) and (iii) we described that the variation of hpF2 and foF2, 
including the seasonal anomaly, cannot be accounted for only, by the geomagnetic 
effet unless other geographical or solar factors are taken into consideration. It 
is the next problem to study what this factor is. 

(iv) That the seasonal variation of electron density may be explained by 
both the vertical electron drift and the variation of the temperature, was suggested 
by Wiss [21]. We briefly touch this. It is considered that the temperature variation 
acceptable physically is very small [5]. If the temperature is expressed by (12) the 
production of the electron is given by 

q= An dS.( = )O)exp| An Tod'seex( a) * | 


0 


l= a , (16) 
where A, S..,m, g, and k represent respectively the absorption coefficient of the par 
ticle, the total solar flux at the top of the atmosphere, the mass of the particle, the 
gravitational constant and Boltzmann’s constant. #) and T) represent respectiveley 
the particle number and the temperature at the reference level. The maximum 


production of electron (q,,) is given by 


dn = AndS.exp| — ( 14% in: (17) 


with An,T,secz=1+r7, (18) 
where 2, and 7, represent respectively the number density of the particle and 
the temperature at the level of g,,. The ratio of gq, in summer to that in winter 
is given by 

(dm)s __ COS (2 +r)e teh Ti 

(Gon ao Cos yu(l-br),e vt ie pe (19) 
where suffixes refer to the summer and winter. The value of Z is 33.1 x 10-°K/cm 
16.5x10-°K/em and 188° x10-°K/em_ for nitrogen molectule, nitrogen atom and 
oxygen atom. Since y may be much less than (probably 3 x 10-°K/cm-5 x 10-°K/em) 
there still be (q,,);>(dm), by a factor of 1.5 or 2.0 in middle latitudes even if in 
summer 7 becomes twice that in winter. Hence the decrease of 4m resulting from 
the increase of temparature is small. 


Next the daily variation of the electron density due to variation of the tem- 
perature without the electron drift, is given by 
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ON Mi OT. -( 1 ae (20) 


a ee ape Ve te 5 
We assume the daytime departure of temperature from the mean (the value at 
06500™) is 10%, then T in the bracket in right side of (20) may be taken to be 
constant. If the temperature varies as follows 


then ar = T,o-cosot. 
Since 7,/7T,=1/10, a ae ~7.3X10-°/s at 06400™ and 7 oh =0 at noon. Hence 


ror", 
rT of 3S less than #. Thus the solution of (20) is not much different from that 


of the equation without the tempererature term. In consequence the temperature 
variation physically acceptable gives little influences on the electron density. 

(v) It appears that the difficulty in (ii) and the phenomena in (iii) in this 
section cannot be accounted for by combining eftect of both vertical electron drift 
and the temperature variation. In order to explain these matters simultaneously, 
the following is suggested. That is, there must be a seasonal variation of the height 
distribution of the particle responsible for the F2 region. For example, we pre- 
sume that various kinds of particle distribute with the common scale height, in 
other words, they are in a mixing state, when the sun is overhead, (for large 
cos zx) and on the other hand, they are probably in a diffusive equilibrium state with 
different scale height one another for small cosy. The distribution variation of 
the particle is due to the deeper penetration of the solar radiation into the 
atmosphere for large cosz. In winter the particle distributes in a diffusive equi- 
librium in high latitude and in summer the partcle approaches to the mixing 
state. Thus the variations of the distribution is larger for higher latitude. In low 
latitude, including the equatorial zone, the sun is always overhead and hence par- 
ticles are always in a mixing state, though the state may be not the same as that 
in summer in middle latitude. The distribution variation in the F2 region takes 
place also when the original variation of distribution is present below the F2 region 
(e.g. near the E region). If the #2 region is formed by ionization of the particle 
which increases its number density at certain heights when the distribution varies, 
the phenomena mentioned in (ii) and (iii) are easily accounted for. The reason is 
as follows. If the particle distributes in the atmosphere where the scale height in 
creases linearly with the height gradient «, the maximum electron production is 
representedd by 

Im= An,0S..exp[—(1+a)], (22) 
with An,H,sec7=1+a. (23) 
In certain latitudes we suppose that S., 7, 4, x, A are all constant. Then the 
height of q,, depends on the product 2H, and its magnitude depends on 2,, the 
number density at the level of the maximum production. If the number density 
of the particle at certain level incraeses which means the increase of scale height 
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(and a), when the atmospheric state changes from the diftusive equilibrium to the 
mixing, the level of the maximum production of the electron rises ((iii) in § 8) and 
the number of the electron production is reduced due to the decrease of m, accom- 
panying increase of H, as seen in (23). Since HM; depends on the mass of the 
particle and the temperature, a fairly large reduction of g, from that in the nor- 
mal diffusive state, about by a factor of 2-4, is expected in summer. So that the 
anomalous seasonal variation of f,F2 may be accounted for. 

Now if the elevation of the maximum production is possible owing to the 
variation of the distribution of the particle, (d) in (ii) is easily accounted for, 
because in winter the distribution approaches the diffusive equilibrium state more 
and more as the latitude increases and the number density at the same height 
becomes smaller. Since the electron production decreases much more in lower 
latitude, compared with that in a diffusive epuilibrium, electron density decreases 
and deviates from the line of fpF 2c cos? z. Thus (e) is explained. If this variation 
of distribution also present in the F1 region also (f) is possible. If during s.s. 
max. the change of the distribution to the mxing state is more promoted than dur- 
ing s.s. min. both reduction of the electron production at maximum level and the 
elevation of its level will be amplified. Thus (a), (b) and (e), are explained. In 
the equatorial zone below 20° there is little change of the distribution of the par- 
ticle so that variation of seasonal electron density obeys to such a way as that 
of the E and F1 region. This is consistent with the observed data. 

In consequence, it appears that there are two main causes for the variation 
(including so-called anomalies) of f,F2 and hyF2. One is the geomagnetic distortion 
which is brought by the vertical electron drift and the other the geographic effect 
which is here considered as variation of the distribution of the particle responsible 
for the F2 region. The former is predominant near the magnetic equator, and 
the latter, in middle and high latitudes. Hence near the equator, the electron 
density is decreased even by a factor of 4-6 during s.s. max., compared with that in 


the case of no mixing and no drift. (see Fig. 2) 
8. Discussions 


(i) It is a very important question whether the solar diurnal and semi- 
diurnal drift are present in the F2 region. The lunar and solar semi-diurnal drifts 
have already been discussed by Martyn. Hirono [22] also calculated the lunar 
semi-diurnal drift, taking the anisotropic conductivity of the ionosphere into con- 
sideration. According to him, the phase of the drift is constant from pole to 
equator over two hemispheres, and he predicted that in the solar semi-diurnal drift 
the relationship of the phase with latitude is closely similar to that of lunar case. 
This result coinsides with that concerning the phase shown in § 3. As to the solar 
diurnal drift the cause has not yet been known. But the presumption from the 
fact that there are the diurnal and semi-dirnal components in the S; variation of 
the earth’s magnetic field, in which the former is larger than the latter, leads to the 
presence of an electric current system responsible for the diurnal component. It 
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is considered that there may be an electric current (or electric field) in the F2 re- 
gion, which is associated with the above current system, and such a observed diurnal 
drift is set up. If the electric field which is estimated from the S, variation is 
present in the F2 region in magnetic equator, upward drift in daytime throughout 
the year is set up. [20]. 

As one of the casuses of the drift, it is expected that the electron moves along 
the magnetic line of force of the earth owing to the diunal wind which varies with 
season. For this purpose, the wind which blows from north to south in summer 
and opposite in winter in the northern hemisphere is necessary. However, Briggs 
and Spencer [23] indicated that many observations in the world show the reverse 
winds respectively in summer and winter to that which we need. Hence the direct 
drift of the electron due to the wind is not consistent with the present observed data. 

(ii) The latitudinal distribution of f,.F2 in Dec. shows that it obeys the law 
of f,F 2<cos?z. We inferred from this fact that the decay law of the electron is 
of an attachment type. Morever in § 7 it is found that decay of the electron by 
the attachment law is appropriate to account for the daily variation of f)F2. 
Besides these, there is another important phenomenon in favour of the attachment 
law ; i.e., the rate of the increase of electron density in the F2 region from s.s. max. to 
s.s. min. differs from that of the Flregion. This difference is removed if the decay 
law in the F2 region is of the attachment type while that in the F1 region is of recom- 
bination type. Further evidences to adopt the attchment law were shown by Yone- 
zawa [24]. So that it seems that the attachment law should be applied to the F2 
region. 

We can presume the magnitude of # in the F2 region in winter in the follow- 
ing ways. . 

(a) The variation of h’F2 in winter shows the conspicuous semi-diurnal 
variation. This implies that variation of the height is not at all due to the penet- 
ration of the solar radiation only, but the semi-diurnal drift influences on the height 
variation. In order to permit the effect of the drift 6 should be less than 10-‘/s. 

(b) The time lag of the maximum density of the day from the noon is about 
3 hours for 2.7x10-*/s. To get the observed time lag # will be 5x10-°/s-7x 10-°/s. 

(c) The ratio of the electron density at sunrise to the maximum of the day 
is about 1/10-1/15 (observed value). According to our calculation this ratio is 
obtained if 8=3x10-°/s-5x10-*/s. 

(d) The magnitude of the electron production is about 80/cc/s for 2.7x 10-°/s 
when the maximum electron denity is taken as 1.5x10°/cc. If the production has 
the value near to this, 8=3x10-°/s-5x10-*/s is hopeful. 

From these points of view, the value of fin the F2 region in winter may be 
3x 10-3/s-5 x 10-*/s and the height corresponding to this value may be 250km-300 km. 

(iii) We suggested a seasonal variation of the number density of the particle 
at any height in the F2 region. Nicolet and Mange [25] indicated that there may 
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be a downward movement of particles in the F2 region, in relation with the down- 
ward movement of the oxygen atom in the lower ionosphere near the E region. 
The diffusive separation may take place at about 160km level ({25], [26], [27]). Hence 
there may be possible that the particle which contributes to the F2 region, distri- 
buted in diffusive equilibrium in winter, changes its distribution owing to the 
mixing in summer and increases the number density at a certain level of the 
electron production. 

Above the height of 200km the main particles are probably O, Ne and UN. 
For example we assume that the F2 region is attributed to the ionization of Nz 
and this particle distributes in a diffusive equilibrium state from the height of 
160km in winter. Let the temperature at this height be 430°K and above that 
height the temperature increases linearly with the gradient of 3°K/km, then in 
summer, owing to the transition to the mixing state, O and N, which is considered 
to be present above 200km [28], move downwards and in turn, N. upwards. Thus 
the number density of N. may increases by a factor of 5-8 at 250km. So that 
the level of maximum production rises higher than one scale height. If the 
particle responsible for the F2 region is O, this phenomenon is reversed. In the 
case of N, its variation of number density is very complicated because the dis- 
sociation of N. must be taken into consideration. 
9. Conclusion 

It is found that the anomalous variations in the F2 region are ascribed 
both to geomagnetic and to geographic causes. Of anomalies, the trough in mag- 
netic equator can be satisfactorilly accounted for by the diurnal electron drift. The 
reason for geographic control of /,F2 in writer is also described. Other variations, 
i.e. seasonal variations of f.F2 and several behaviours of h,F2 and f,F2 pointed 
out newly are explained if we assume that the height distribution of particle 
varies seasonaly. 
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Geomagnetic Disturbances* 
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Abstract 


The existence of some correlation between the geomagnetic 
storm and the ionospheric F2 layer disturbance has long been in- 
vestigated since the early stages of ionospheric researches. Recenty, 
this problem has been studied by many investigators not only from 
geophysical interest but also under the necessity of improving radio 
communications. 

The main course of investigation on F2 disturbance, which 
has been studied by the author for the last three years, is sum- 
marized as follows: 

(1) The separation of variations in /’F2 and h'F2 associated 
with geomagnetic storms into two parts depending on the universal 
time and local time, the so-called Dsi(F2) and S,(F2), respectively, 

was clarified. 

(2) The devlopment process of the F2 disturbance was in- 
vestigated in respect to each of D,,(F2) and S,(F2). It is concluded 
that S)(f’F 2) exists in the earlier stage of geomagnetic storm and 
D,(f’F2) grows in the comparatively later stage. 

(3) The latitudinal and seasonal characteristics were shown. 

(4) Investigations were made in the behaviour of the F2 
layer according to the different levels of geomagnetic activities. 
Results showed that the F2 layer disturbance was generally in 
proportion to the geomagnetic activities, except for D,,(f/’F2) vari- 
ation in winter. 

(5) As to the successive yearly variation of. ionospheric dis- 
turbance, it was found that the phase of S,(F 2) would lag behind 
according to the decrease of sunspot activity. 

1. Introduction 

The existence of irregular variations in the F 2-layer associated with geo- 
magnetic disturbances had been investigated in the ealier stages of ionspheric 
researches by E. V. Appleton, etc. [1] and L. V. Berkner, etc. [2] By gradual steps, 


* Previous reports (3) (4) and (5) on investigations by the author are collected in this paper, 
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the morphology of ionosphere disturbance was being clarified, but these investi- 
gations were intermitted by World War II until some four years ago. Recently, this 
problem has been studied by many investigators [3] from both geophysical interest 
and practical necessity. 

The first step to do was the statistical separaiton of variations in f’F2 and 
h'F2 into parts depending on the storm and local time (longitude), respectivly, as 
was done in the case of variations in geomagnetic forces during geomagnetic 
disturbances. 

After this was done, the morphology of F2 disturbance was clarified without 
difficulty and some foundation was laid for the warning of radiocommunication 
disturbance. 

2. Separation of F2 Disturbance 
into Dy(F2) and Sp(F2) 

The success of separation into 


Washington Summer (1944~1951!) 


D,, and Sp parts of the geomagnetic 
variation suggests the same -suc- 
cess in the case of F2 disturbance 
because of close interrelationship 
existing between them. Under this 
consideration, we are applying the 


same method to the ionospheric < 


data in this paper. Pithdor 
As the ionospheric quantities Vaca EN NOY AACA CaO 
there are taken the values of 
0 
pate (where <4f°F2 represents Fig. 1 Ds; (left) and Sp, (right) parts of the F2 
Laas bts vs disturbance at Washington and Wak- 
the deviation of f°F2 from the AfoF2 


kanai. Variations in by dotts 


AfoF2 
monthly median value, and /°F2* FoF 2 * 
; and Ah’F2 by lines. 
represents the monthly median 
value) and 4h'F2 for many stations in different seasons. 
As shown in Fig. 1 quoted from my earlier report [4], we can easily find that 


the F2 variation contains two variable parts like D,, and Sp, during the magnetic 


storm. 
3. Characteristics of F2 Disturbance 


Summer 3. 1. Developing process of F2 disturbance 
The development of the D,,(F'2) part is shown 
= swormTine in Fig. 1. As shown in Fig. 3, D,, disturbance of 
Winter h'F2 begins gradually without delay after SC in 


spite of the fact that /’F2 variation has not an 
appreciable change immediately, and then large 

ie Zi earings Pe Oy me variation follows 8~10 hours after SC. This fact 
Washington (1944~51). shows that the ionospheric disturbance is not of 


24 StormTime 
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propagation phenomenon from high latitude toward low latitude, but it begins 


simultaneously with SC. 
The development of the 


S,(F2) part is shown, under the 
approximation of the sine curve, 
in Fig. 2 at Washington (1944~ 
51) as described in the previous 
report [5]. By combining Fig. 1 
and Fig. 2 we can reproduce the 
original disturbances which 
have begun with SC occurring 
at different local times as 
shown in Fig. 3. 

On the other hand, in 
order to study the development 
of S,(F2), without assuming the 


Washington Summer (Dstt So) and MeanValucs 


0 Loca! Time 
of, 4 Slorms 


% Washington Wiriter (Dst+Sp) and Mean Values 
20 a 
10 Local Time ‘ 
° 
teanot 9 Storms mean 4 | 
4 


oes ESS 


i) 
+ whe 


constant phase, we can use the 
storm time diagram of /°F2. 


Fig. 3 Development of F2 disturbances in f°F2 at 
Washington (1944~51), having their sudden com- This was done for Washington 


mencement at different local times. f° F2 [3] in summer and winter, 


as shown in Fig. 4. 


We can see that D,,(F2) begins with a large change at about 8~10 hours and 
the maximum amplitude of S, occurs 6~12 hours after SC. Then it is found that 
S> is large in the early stage of the magnetic storm, but D,, predominates later on. 
It seems that the latter has not so good flexibility as the former according to the 
magnetic variation. As seen in Fig. 4, the phase lag of Sp is clear especially at 
the last phase of the disturbance. This fact was also pointed out by Mr. T. Oba- 
yashi [3] in his study of individual behaviour of F2 disturbance. 

The above development processs of F2 
disturbance will be helpful to the prediction 
of radio communiction disturabnce. For in- 


Summer 


stance, Fig. 3 may serve as an important key 
to forecast the world-wide distribution and 
development of disturbance in regions in 
which is measured the lower field intensity 
of radio waves. 


Storm Time 


3.2. Latitudinal and seasonel character- °0 6a 


iS 2 
istics of F2 disturban cate Local Ties 
en Fig. 4 Storm-local time diagram fo Sp 


a. Latitudinal characteristics (f°F 2) at Washington (1944~51) 
The data used concern /°F2 observed in summer and winter. 


at Fraserburgh, Washington, Wakkanai, Yamagawa, Huancayo, Singapore, Brisbane 
Falkland Island, Hobert and Macquarie Island, during the summers and winters of 
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from 1949 to 1951 (Fig.5). 

The results are expressed in geo- 
magnetic latitude local-time diagrams of 
S,(f°F2) and geomagnetic latitude-Dst(/’F 2) 
diagrams, as are shown in Fig. 6 (a) and 
(b). Here, D,.(f’F2) are the mean values 
of D,.(f°F2) during the period of 48 hours 
subsequent to SC in the stormy time. 

In this figure is shown that the phase 
of S,(f°F2) lags from summer to winter 
hemisphere. It is strange that this phase 
changes gradually to about 180° without an 
equatorial node. D,,(f’F2) illustrates the 
well-known fact of the summer decrease 
and winter increase in /’F2 with the ex- 
ception of the singular area of D,,(/’F2) 
which exist at about 50° geomagnetic lati- 
tudes as shown Fig. 6. 

b. Seasonal characteristics 

As are shown in Figs. 5 and 6, the 
D,.(f°F2) variation has a general tendency 


in summer and increase in 


to decrease 
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Fig. 5 (a), (b).F2 disturbances of f°"F2. When 
the northern hemisphere is in summer 
(a) and in winter (b). Shaded areas 
indicate the negative values. of 


Sp( fF 2). 


Fig. 6 (a), (b). Both (a) and (b) show 
f°F2 disturbances in percent- 
age with respect to geomagne- 
tic latitude and local time, in 
different solstitial seasons, 
respectively, as described in 
Fig. 5. The left show Dst 
part, and the right S, part, 
both being thea verages of 48 
hours in stormy time. 
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Fig. 7 The average of 24 hourly mean values of 
A 0 


fu Fo on the five international monthly 
geomagnetically disturbed days and quiet days 
at Washington and Kokubunji. Thick lines 
are for disturbed days and thin lines for qniet 


days. 
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Fig. 9 (a), (b). 


winter. 
Fig. 7 shows the average of 
AZoE: 


24 hourly mean values of fv FQ% 
on the five international monthy 
geomagnetically disturbed days 
and quiet days both at Washing- 
ton and Kokubunji. From this 
figure, we can see the seasonal 
and yearly variation of D,.(/°F 2), 
but the latter does not well 
appear, particularly at Washing- 
ton, owing to the values divided 
by median values. 

The average curves of D,, 
(f°F2) for about 4 years from 
July, 1949, to February, 1954, and 
D,,(h'F2) for about 2 years in the 
same period are shown in Fig. 8 
(a) and (b). The tendency that 
h'F2 is elevated on geomagne- 
tically disturbed days and is de- 
pressed on quiet days is witness- 
ed in common at both stations. 
It is noticeable, however, that the 
curves of disturbance days in f°F2 
have two minima near the equi- 
noxes and the curves of quiet 
days have two maxima, and 
both curves are mirror images 
of each other. These maxi- 
ma and minima are not 
attributed to the equinoctially 
maximum activities of geo- 
magnetism, for the reason 
stated in the following sub- 
section. Hence, it may be con- 
sidered that Fig. 8 expresses 
the seasonal tendency of the 


Deviations of f°F2 and h'F2 at different degrees ionospheric disturbing sensi- 
of activities from the groups of the calmest magnetic bility. 
activity, are summarized in three principal quanti- 


ties: the magnitude of Dst, Sp and the phase of S 
respectively, 


Dy 
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3.3. On relations between the grade of geomagnetic activity and F2 disturbance 

The ionospheric variations are investigated on the basis of the data of 
Washington over the years 1949~1953 divided into three groups of the year. 

The geomagnetic activities are expressed by SKp, which is the summation 
of Kp’s of the whole day. The corresponding F2 values used here are the deviations 
of the mean values of f°F2(h'F2) which belong to the groups of days of 13-<2Kp 
S22<+, 23-<YKp<32+ and SKp>=33-, from the values belonging to the quietest 
group YAp<12+. 

In order to show the principal characteristics of these deviations, noncyclic 
variations [6] are subtracted from the above deviations and they are smoothed by 
means of three hourly running averages. The three main characteristics, that is, 
the mean values over 24 hours, the range of diurnal variations, and the centre time 
between the times of maximum and minimum values of the curves, for each group 
of days in the three seasons, respectively, are summerized in Fig. 9 (a) and (b), and 
an example of the original curves f’F2 in case of )Kp=33- is shown in Fig. 10. 
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ra 6 6 
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v c c 
3 
E : bE 
Go “6 RR 1s 0 6 2 8B i ae 750° 
Local Time Local Time Local Time Sunspol Number Rz Sunspot Number Ra 
Fig. 10 Example of the original curves in Fig. 11 Variation of the phase of S, and the 
f°F2 corresponding to the case of level of Ds, at Washington accord- 
Kp=33-. ing to the sunspot number Ry, in 


different seasons. 


Fig. 9 shows that F2 layer is generally disturbed in proportion to the geo- 
magnetic activity except for the D,,(f’F2) variation in winter. 

At Akita and Brisbane, these former characteristics are also accepted. 

3.4. Yearly characteristics of F2 disturbance 

For using os which is divided by the monthly median valus as shown in 
Figs. 7 and 9, the yearly variation of #2 disturbance is not remarkable in the 
amplitude of seasonal variation as described above, but 4/°F2 is conspicuous by its 
successive yearly variation in proportion to the solar activity. 

The same character exists as well in S,(f°F2), and additionally Sp has a 
tendency to advancement of the phase according to solar activities. (Fig. 11) 
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In summer and equinoctial seasons, the variations of /’F2(h'F2) decrease 
(increase) continuously as geomagnetic activities increase. In winter, however, the 
behaviour of F2 layer is not so simple, for f’F2 rises from slight to moderate 
magnetic activity and h’F2 monotonously increases as 
illustrated in Fig. 12. This behaviour is considered to be 


Washington 


3 1959-50 
the general nature of the disturbance in winter, as was 


discovered by L. V. Berkner, [2] etc. at Wathroo during 
the epoch of the larger sunspot number. 
4. Concluding Remarks 

The above characteristics of the F2 disturbance 
associated with geomagnetic disturbances are recognized, 
without any remarkable contradiction, by other investi- 
gators somewhat different methods and data. There will 


be some time before the complete morphology of F2 fig 12 Diagrams of 
disturbance is accomplished. And yet on this ground will Ds(f°F2) and 
be founded the theory of ionospheric disturbance and the Dud E a) 

prediction of radio disturbance. In closing, the writer wishess to express his 
sincere thanks to Dr. H. Uyeda and Prof. T. Nagata for their continued interest 


and encouragement given him throughout this work. 
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Geomagnetic Distortion of the F2 Region 
on the Magnetic Equator 


By Motokazu HIRONO and Hiroshi MAEDA 
Geophysical Institute, Kyoto University 


Abstract 


The direct relation between the geomagnetic S, variation and 
the vertical electron drift in the F2 region on the magnetic equator 
is examined. 

It is shown that the electric field in the F2 region accompanied 
by S, electric current produces the vertical drift which is sufficient 
to interpret the main features of the anomaly of the F2 region on 
the equator. It is to be noticed that the main term of the drift 
velocity is diurnal. The daily variations of the maximum electron 
density and its height in the F2 region are calculated under con- 
sideration of the vertical electron drift for the reasonable distribu- 
tion of decay coefficient with altitude inferred by observed results. 
The calculated F2 daily variations have a striking resemblance 
with those observed near the magnetic equator. When the ion 
production takes its maximum value at about 200km, there appears 
a lower secondary maximum of electron density which agrees well 
with the observed F1 layer. 

The change of characteristics of daily variations of the F2 
region with the sunspot-cycle is likely to be accounted for by a 
slight shift of the phase of the drift. 

1. Introduction 

Analyses of ionospheric data and theoretical researches by Martyn [1], Mitra 
[2] and Weiss [3] have established the existence of solar semi-diurnal components 
of tidal origin in the maximum electron concentration N,, of the F2 region and 
its height ,,, and further studies have been made by K. Maeda [4]. 

We tried to find a direct relation between the daily variation of the earth’s 
magnetic field and the vertical electron drift in the #2 region on the magnetic 
equator. 

McNish and Gautier [5] have attacked a similar problem, and an interesting 
relation between these data is shown. Their theoretical consideration is, however, 
confined to the vertical drift due to the electric field by the daily variation of the 
magnetic field, whereas the drift due to the electrostatic field is not taken into 
account. According to our results, the latter is so much greater than the former 
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that we disregarded the drift due to electromagnetic induction. 
2. Observed Data of F Region near the Equator 

The world-wide distribution of N,,(F2) which is closely related with the 
earth’s magnetic field was studied by Uyeda and others [6] and by Appleton [7]. 
Further studies have been made by Bailey [8], Uyeda [9] and Aono [10]. According 
to their results, for noon conditions, there is a belt of low values of N,,(F2) circling 
the Earth and centred roughly on the magnetic equator, and the maximum value 
of N,,(F2) is attained roughly at about the 20° north and south geomagnetic 
latitudes. The data analysed in detail came from eight widely distributed stations 
close to the magnetic equator and about the 20° geomagnetic latitudes (which will 
be referred to as Eq. and 20° stations respectively hereafter) for which the 
ionospheric characteristics at sp-maximum (sp=sunspot) and/or sp-minimum are 
well observed. All information was extracted from the “Ionospheric Data” of 
C.R.P.L., Boulder, except that of Palau. The stations used are shown in Table 1. 


Table 1 
Station Geographic Geomagnetic Geomagnetic Times 
i iat: Long. Lat. Long. Dip Standard 

Ibadan 7°.4N 4°.0E 10°.6N 74°.7E —5° 0° 
Tiruchy 10°.8N 78° .8E 1°.0N 148°.4E 4° Local 
Leyte 11°.0N 125°.0E 0°.28 193°.8E 4° 135°E 
Palau I. 7°.3N 134°.5E rages 156°.3W 2° 135° E 
Christmas I. 1°.9N 157°.3W 2° .2N 88°.7W 6° 150° W 
Huancayo 12°:0:S 75° .3W 0°.6S 6°.2W 2° 75°W 
Maui 20°.8N 156°.5W 20°.9N 88°.1W 40° 150°W 
Trinidad 10°.6N 61°.2W 22°.0N 8°.3E 43° 60°W 


The values of N,,(F2) are calculated from f°F2 and h,, from the value of 

(M3000)F2 which will be denoted by M, using a formula by Shimazaki [11], ie., 
h» = —176+-1490/M (km). 

Various daily variation curves are shown in Figs. 1, 2 and 3. To see the change 

of characteristics with the sp-cycle, lack of data is to be met with at times. Hence 

the data of the years as close as possible to sp-max. or sp-min. were used. 

For Eq. stations the daily variations are averaged from one year and for 20° 
stations the N curves averaged from the Northern solstice, May, June, July and 
Aug., the FE curves from the equinoctial months, Mar., Apr., Sept. and Oct., and 
the S curves from the Southern solstice, Nov., Dec., Jan. and Feb. For Huancayo, 
h'F 2, h'F1 and N,,(F1) are illustrated as well, and some data analysed by Martyn 
[1] are shown in addition. 

It can be seen from the figures that the daily variations at Eq. stations bear 
a good resemblance, moreover they are very similar to those observed at summer 
conditions in the middle and high latitudes. The values of N. and hh, ata 


™ 


stations increase with the increasing sp-numbers. ‘The values of N,A(F2) and 
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Fig. 1 Daily variations of Nm(F2) and hm(F2) at stations close to the magnetic equator. 
The values of R denote the mean sunspot number for the year. 
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Fig. 2 Daily variations of Nm(F2) and h»(F2) at stations 
of about 20° geomagnetic latitudes. Dotted lines 
refer to the northern solstice (N), dashed lines to the 
southern solstice (S), and full lines to the quinox (E). 
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Fig. 3 Daily variations of Nm(F2), Am(F2), h'(F2), Nn(F1) 
and f'(F1) at Huancayo in years of different sun- 
spot numbers. 


N,,(F1) at noon of sp-max. at 
Huancayo are about 2.4 and 1.6 
times of those at sp-min. The 
forenoon maximum of W,,(F 2) 
is smaller than the afternoon 


‘maximum of WN,,(F2) for the 


sp-minimum period. The for- 
mer increases in relation to the 
latter with the increase of sp- 
number. At Huancayo, especial- 
ly, the former is greater than 
the latter at sp-max. 

The daily variations of 
h,,(F2) at Huancayo and Maui 
were subjected to a harmonic 
analysis and the amplitudes P,, 
and phases ¢, when maximum 
value was attained for mth 
term are shown in Table 2. 

It can be seen from this 
table that at Huancayo the 
diurnal term is much greater 
than the semi-diurnal one, 
whereas at Maui the relation 
is reversed. 


Table 2 
@ — P 

: Bis Pr t Pe, VoVG P. t 

St t Y 0 1 L 2 2 3 3 
ation ear Season (km) (km) | (hr.) | (km) (hr.) | (km) | (hr.) 
Huancayo 1942-44 mean 350 61.9 14.2 Ia.c 10.7 14.6 ace 
1940 mean 394 87.8 14.9 Ley, 9.7 16.0 ee 
1948 mean 429 110.6 15.7 16.8 9.4 8.0 pss 
Maui 1945 N 327 28.0 12.0 SUE 0.6 4.6 4.6 
¥ E 306 9.2 9.6 2157 10.9 16.8 Soo 
hs % S 294 14.2 6.3 8.6 10.7 25.6 5.6 


3. Vertical Electron Drift of the F Region on the Magnetic Equator 
(i) The motion of charged particles and the electrical conductivity of the 


ionosphere have been discussed by Hirono [12], [13], [14] (hereafter referred to as 
I, Il and III respectively). 


When the electric field attains such a value that the vertical electric current 
vanishes, the velocity of the vertical drift, which is the same for ions and electrons 


and measured positive upwards, is given on the magnetic equator, according to 


Ill, by 
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W=f(h)E,/H, (1) 
where H is the magnitude of the main magnetic field, and FE, is the eastward 
component of the electric field, and f(h)~p;/(p;+p,) with p,=(1/m,){»,/ve2+o,7)}, 
(y=i,e), where v; and v, are roughly the collision frequencies for ions and electrons 5 
, and w, are their frequencies of spiralling round the lines of magnetic force ; m; 
and m, are their masses. For this expression, the particle concentration of nega- 
tive ions is taken as effectively zero according to Bates and Massey [15]. 
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Fig. 4 Atmospheric model and distributions Fig. 5 Variation of f(h)= W/(E,/H)=p:/(pit+ pe) 
of electron density with altitude. with altitude. 

Curve A denots the Chapman layer, 

and curve B is inferred from rocket 

observations [21]. 


We adopt an atmospheric model inferred by rocket observation [16] as shown 
in Fig. 4 to which the following calculations refer. As illustrated in Fig. 5, f(h) is 
nearly unity above the height of about 120km and decreases rapidly in the lower 
region. When the vertical current is not much greater than the horizontal one, 
Equation (1) holds approximately for the vertical electron drift in the F region. 

(ii) According to the dynamo theory, the electromotive force of the S, cur- 
rent on the magnetic equator is the electrostatic field. It is shown in II from the 
theory of potential that the eastward component £, of this field is almost constant 
between the £ and F regions. 

The electrical conductivity of the ionosphere in the E-W direction is expres- 
sed, as shown in I, by 

63=0;(1+027/0;"), 


a nite Vv, nN; Vv; 
where a,=e( Pat Sai teh ened Ota oF hice) 
€ 4 Vet; 
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where , and , are particle concentrations of electrons and ions, and e denotes 
the absolute value of the charge of an electron. Let the eastward electric currents 
integrated with height in the E and F regions be I, and i respectively, then we 
have approximately 
h=K,:E,, L=kK,-Ey,, (2) 
where 
K,=\ o,dh, Ks = | o3dh. 
Ez F 


According to II, III and H. Maeda [17], Ki>Ke; therefore it follows that 
L>h. 
Put I,+h=-J, and let the magnetic variation due to J at the lower edge of the E 
region be 0H, then we have 
6H=2rI~2rl,. (3) 

Considering the fact that about 80 per cent of the daily magnetic variation 4H on 
the ground is caused by J and that this magnetic field due to J increases upwards 
and the amount of the increase will be about 20 per cent at the lower edge of the 
E region [18], so that we have approximately 

4H~0H~2rK,- Ey. (4) 
When the values of 4H and K, are known, we can determine the value of EF, using 
this relation. 

(iii) M. Hasegawa and H. Maeda [19] have estimated the daily variations of 
electrical conductivity of the S, layers (mainly the E layer), and as a result of this 
analysis the daily variations of electric field and the zero level (C)) of the S, varia- 
tions have also been estimated. Further, H. Maeda [20] has made his calculations 
on Huancayo’s geomagnetic data for a period of 12 years, and the results for the 
mean states of sp-min. years, 1922, 23, 32 and 33, and sp-max. years, 1926-29, are 


as follows: 
4H=Cy+ dC, sin( 2d + «,) (10-°e.m.u.). (5) 
~ Average pe hl Fi ek | fod | onl Sais a 
Years pe per Co | Cy | ey | Cy | Eo | Cs | &3 
1922, 23 z as ; 
32, 33, 8 (min) | 15.1 | 32.5 | 280°.2 | 15.8 | 108°.4] 5.8 | 306°.0 
1926, 27 : “a ae : 
“98, 29| 71 (max) | 24.3 | 52.4 | 274°.1 | 25.7 | 103°.9| 10.9 | 302°.2 
K,=K) {1+ d)7,,cos(zd+6,,)} (e.m.u.), (6) 
v1 | 6, | Y2 | bs 


1.15 | 180° | 0.42 | 0° 


(almost independent of sp-number) 
where A is the local time is angular measure from midnight and 4H is measured 
positive northwards. The daily variation of K,/K, is shown in Fig. 6. 
On the magnetic equator, Hall conductivities d,y and oy, vanish and Ohm’s 
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law holds in its original form, so that x, 
the results estimated by the method 24 
of the above authors are theoretically 20 
adequate. The noon value of inte- 
grated conductivity of the E layer cal- 
culated as a simple Chapman layer of 
scale height=10km, maximum electron 
density at noon=1.5x10°/c.c. (for sp- 
min.), height of maximum ion production 00 04 08 2 16 
=100km as shown by curve A in Fig. Fig. 6 Daily aise paar conductivity 
4, is 3.4x10-'e.m.u., so that from the in units Ko. 
relation (6) K, becomes 1.3x10-"e.m.u. For the sp-max. period mentioned above we 
take K,=1.65x10-'e.m.u. corresponding to the maximum electron density of 1.9 
x10°/c.c. for the E layer. 

_ The electrostatic field E, expressed in Fourier series does not include a 


20 00 


constant term, since this is derived from the potential. Therefore we have the 
following expression of FE, and Table 3 in connection with relations (4), (5) and (6), 
using the value of K, estimated above: 


E,=y:sin(A+9,)+yesin (24+62) (e.m.u.). 7a) 
Table. 3 If we take the value of H=0.27 
auss at the height of the F 
Years / yy | 61 | 2 | 82 & , & 
| region over Huancayo, from (1) 
sp-min. 3.3410? 280 0.3610? 176 tad (7) welliave: the dollowtir 


0.8110? | 17h" 


sp-max. 4.19X1F ) 265° utr eel iis: 
W=f(h){W,sin(A+9,)+ Wesin(22+-9:)}, (8) 
where 
Sr Adan neh iosiepdateperitd: 
W,.=0.13 x 10°cm/sec. (9) 
W1=1.55 x 10°cm/ sie for sp-max. period. 
W.=0.30 x 10°cm/sec. 

The daily variation of W for sp-min. 
eit at an altitude of 300km is shown in 
= : Fig. 7. The value of K, is influenced 
» 
3° to some extent by the assumed elect-. 
ss ron density distribution with altitude. 
a If the electron density above 100km 

00 04 08 12 16 20 0 is 1.5x10°/c.c. and below that level the 


pee er yer : same as the above mentiond Chapman 
Fig. 7 Daily variation of vertical electron drift 


velocity W for sp-min. at a height of 
300km referred to curve A in Fig. 4. cent. This variation is very small. 


On the other hand, when we adopt an electron density distribution with 
altitude inferred by rocket observation [21] as shown by curve B in Fig. 4, A, takes 


layer, K, is increased by about 10 per 
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one second of the above mentioned values, i.e., 6.5x10-%e.m.u. for sp-min. and 
8.3 107Se.m.u. for sp-max. If these values of Ky are used, the electric field E, and 
drift veloity W become twice as great as those represented by (8) and (9). 

Effective conductivity o; is obtained under the assumption that the electric 
and magnetic fields are completely uniform, therefore if they are not completely 
so, the effective conductivity in the E-W direction may be slightly less than o3. 
Then the decrease in K, will give slightly greater values for W in any case. 
4. Solution of the Continuity Equations 

Let unit of height ® be 50km and z=(h—,)/H, and let z=0 at an altitude of 
300km. In the F region, the local scale height of the atmosphere (F) is taken to 
be the same as the height unit $, but in the E region the local scale height is taken 
to be a function of altitude which may be deduced from Fig. 4. 

The equation of motion for an electron is dh/dt=W which may be trans- 


formed to 
F=137x 10. (10) 
The solution of this equation for a prescribed value of W may be written 
2=2(2, a), (11) 


with z, fixed by initial candi- 
tions at A=Ap. 

In Fig. 8, several curves 
emerging from (11) derived 
from (8) and (9) for the sp-min. 
period are illustrated for some 
values of z. According to this 
figure, it was found that a part 
of the electrons in the F2 
region descends at night and 
enters into the E region. In a 
later section, however, the 
results of calculation show that 
the electron densitiy in the cell, 


00 04 08 12 16 20 00 


LOCAL TIME (HOURS) which descends below 300km, 
Fig. 8 Solutions of the subsidiary equation o =W rapidly decreases because of a 
for some initial values of h. great recombination cofficient 


in the lower region, tending to be several times 10'/c.c. and the electron density 
of the F region does not seem to be appreciably influenced, Thus, #,,(F2) does not 
descend much below the 300km level. 

(i) Theoretical discussion of the effect of the vertical electron drift upon 
an ionized region, necessitates a consideration of the equation of continuity 


ON 2.0 
Tere =J— alN* — on N: W), (12) 
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where WN stands for electron density per c.c., J for ion production per c.c. per sec., 
and « for the recombination coefficient. This equation is to be solved subject to 
the subsidiary condition (11) imposed by the drift. 

In the F region W is estimated to be almost independent of altitude, so that 
we have 

ON. =[I],—[o].N?, (13) 

where d/dit stands for a derivative “following the motion.” In this equation J and 
a are to be evaluated at the height z which is obtained from (11) as a function of 
z and ¢. It is shown by Bates and Massey [22] that a=8x10-"'c.c./sec. in the F2 
region and a=4x10-%c.c./sec. in the F1 region (independent of pressure) in the 
daytime. For the F2 region at night, Ratcliffe [23] obtained a recombination coef- 
ficient of the same order of magnitude as the above daytime value, in the absence 
of the vertical drift influence. Therefore, we take a=a,§(z), where &(z) is taken to 
be unity for z=0 and 50 for z<—1 as illustrated in Fig. 9. For the electrons 
which have descended into the E region 
at night, a third term including the 


6.0 600 
4.0 500 vertical gradient of W appears on the 
— 400 _ right hand side of (13). In the follow- 
Ete A's ing calculation this term is neglected. 
nq“ zo, According to the study by Weiss [3] 
ini : : ae. on the variation of electron density 
a 8 260 — . . . e . 
ae a0) at night, this neglection will be justi- 
-20 200 ~— fied when we confine our attention to 
eg! the daily variation of electron con- 

100 
centration in the F region. For the 

| » ae 5 10 20 30 50 z . . : . 
az) He) distribution of ion production, the 
Fig. 9 Variation of Kan and ¢(zj= 


following two types (a) and (b) are 
ies manana considered : 
(a) The maximum of ion production at noon is located at 300km. 
According to S. Chapman [24] 
I=]F(Z, x), 

where F(z, ~)=exp{1—z—exp(—z)Ch(R+z, z)}, and Ch(k+z, x) is a complicated 
function of R and the solar zenith angle x, and R=150 is used here. After some 
transformation of (13) and by hae N,(a)=(1o/a)'/* and v=N/No, we get 


Sa a =(Fy-lel*, (14) 


where 1/o,=1.37 X 10'(1ya)!”. 

(b) The maximum of ion production at noon zs located at 200km. 

Bradbury [25], Mohler [26] and Bates [27] suggested that the Fl and F2 
layers are produced by the same atoms or molecules ionized by the same ultra- 
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violet, and the two layers are separated by a rapid change of the recombination 
coefficient with altitude. To examine the effect of the drift upon an ionized region 
under such conditions, we take J=J,F(z+2, x). Putting M(6)={hF (2, 9)/a}”, 
v=N/N, and F’=F(z+2, x)/F(2, 0) and after some transformation of (13), we get 


oP = LF, (15) 


where 1/0, =1.37 x 10*{1,F (2, 0)a}1/?. 

(ii) According to the theory of dissociative recombination by Bates and 
Massey [28], a recombination of electrons with ions might take place in the same 
form as that under the attachment law. Though this theory is not yet established 
quantitatively, we shall use the following equation in place of (13) 


dN 
ON =[T1.—I61.N. (16) 


We take @={,)¢’, and the variation of é’ with altiude for z>—1 is taken as shown 
in Fig. 9, and 6,=8x10-*/sec. is adopted from the relation 8; =No%. Ion production 
is taken as type (b), ie. .=1,F(z+2, x). Putting F’=F(z+2, x)/F(2, 0), Nolo) 
=J,-F(2, 0)/B) and »=N/No,we have 


oO = LF —[e" Ly, (17) 


where 1/¢,=1.37x10'8. For the altitudes z<—1, (15) is used instead. 
5. Distortion of the F Region by Drift subject to S, Current 

To evaluate the effect of S, drift on the F region on the equator, we solved 
Equations (14), (15) and (17) following the motion subject to the subsidiary condition 
(11) using the method of Millington [29]. Numerical integration is made at an 
interval of 10° of local time 4, but between the height —1<z<0 the interval was 
taken to be 5°. Using these solutions, the maximum electron density v,, in units 
of N, and its height h,,(=z,,-+300km) are illustrated in Fig. 10. The details of 
the cases calculated are summarized in Table 4. 


Table 4 
Case | Continuity Equ. used DNR RE eng rs a Curve in Fig. 10 
W,(cm/sec) 0, We (cm/sec) Qo slunteatiats 

a (14) with og=1 1.24108 | 280° (115 .3t) | 0.13108 176° (a),Curve A 
2 ne 2.48 ,, - 20. ne a B 
3 (15) with o,=1 Lae 0.13 ,, : (b),Curve A 
4 “ Beko ss "S 026";, > % B 
5 (14) with og=1 lta. 0° ( 6%.0f) | 0 (c),Curve A 
6 Pi: ape ye 180° (185 .0t) | 0 PS B 
7 (17) with o,=1 7 > 280° 0.13 x 10% 2 leiliyg (d) 
8 (14) with og=1/2 io eae a te ae o (e) 
9 (15) with o,=1/1.5 | 2.48 ,, “9 WZ se as (f£) 

10 (15) with o,=1/1.5| 3.10 ,, 265° (125 .3t) | 0 (g) 

11 (15) with o,=1 2.48 3, 250° (13 .3t) | 0 (h) 

12 (15) with o,=1 2.48) 33 220° (16334) 1-0 (i) 
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Figs. 10 (a) to 10 (i) Calculated daily varia- 
: N»(F 2) 
tions of Von(F°2) = N, . hn(F 2), 


0 
Nm 
Vn(F'1) — vm and hn(F1) for 
0 
static (dotted lines) and moving (full 
lines) regions. Every curve is ob- 
tained under the conditions shown 
in Table 4. 


00 04 08 12 18 20 00 
LOCAL TIME (HOURS) 


(i) In the F2 region it is estimated that o,, o, and o, in (14), (15) and (17) 
are close to unity. Therefore Cases 1 to 4 in Table 4 are calculated. In any 
case the corresponding values »,, and h»,, for the static Chapman layer are shown 
in the same figure but by dotted lines. 

The variations of v,, and h,, of these four cases are very similar to those 
observed at Huancayo and other Eq. stations, at the period of minimum and 
moderate sp-numbers. In general v,, is smaller than v»,,, and at noon the decrease 
dv,, from Y,, is between 40 and 60 per cent, and the maximum of »,, is attained in 
the afternoon. The daily variations of h,, are subjected to harmonic analysis, and 
the amplitudes P, and phases ¢, of the terms are shown in Table 5, 
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Comparing these values with those Table. 5 
calculated from the observed data | | 

Case Py ty Py ty Ps ts 
at Huancayo (1942-44) as shown in (km) | (km) | (hr.) | (km) | (hr.) | (km) | (hr.) 
section 2, it is to be noticed that 
Cases 1 and 4 fit well with those 
observed. However, so far as »,, is 
concerned, Cases 3 and 4 are in | 


338) | 62.5:] 13.3 | 7.9) 921°) 13.42) 1.8 
347 | 82.4 | 12.4] 15.9 | 10.2 | 8.3) 2.3 
Oe aliaeG) ee SoSileOrast. Oneal al 
S267} BLS 12200 L067 PL 2.0726 


er wON FE 
oo 
bt 
rar 


better agreement with those observed, moreover in the two cases there appears 
a lower secondary maximum of electron density which may be termed as »,,(F1). 

The daily variations of v,,(F1) and its height h,,(F1) does not show appreci- 
able departures from static Chapman norm because of a great recombination 
coefficient and ion production in the altitudes. At noon we have »,,(F1)N)(b)=2.5 
x10°/c.c., if we take N,(b)=10°/c.c., and at the same time k,,(F1)~200km. Thus, 
the lower layer is very similar to the observed F1 region. 

(ii) The question may arise, whether or not the drift W with other phases 
than those derived from (7) have similar affects upon the ionized region. To 
answer this question Cases 5 and 6 have been calculated. The variations of h,, in 
these cases do not resemble those observed near the magnetic equator. The varia- 
tions of v,, do not show appreciable departures from the static Chapman norm »),,, 
and especially in case 6 we see that v,,~»,. If the value of W appropriately 
tekes less value in case 6, v,, will be greater than »,, by about ten per cent, since in 
this case the motion of a cell is in phase throughout the day with the motion of 
the level of maximum ion production. 

(iii) Case 7 was calculated to find the daily variation under the attachment 
law. In this case the variation of v,, in the daytime resembles those observed, but 
at night v,, is too small compared with those observed. Better results may be 
obtained for different values of 8, but further discussion is not included in the 
present paper. 

(iv) Cases 8 and 9 were calculated to evaluate the effect of the decrease of 
d, and a,. In Case 8, the effect of the drift upon the F2 region is a little decreased 
compared with Case 1. The calculated value of N,,(¥ 2)=v,,+N.(a) at noon for this 
case increases by a factor of 3 compared with that for Case 1, while No(a) by a 
factor of 2. 

In Case 9, the calculated W,,(F2) and WN,,(F1) at noon increase by factors 1.8 
and 1.5 respectively compared with those of Case 4, while N.(b) by a factor of 1.5. 

Though at present we have no firm ground to assert that o,/and-4, are 
closer to unity at the period of sp-minimum, the above mentioned relative varia- 
tions will throw some light on the variation of the F2 region due to the sp-num- 
bers as mentioned in section 2. 

(v) As shown in Table 3 the phase 9, of the diurnal term of £, for the 


sp-max. period is 265° instead of 280° for sp-min. The range of errors for the 
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phases may be several degrees at present stage of calculation and much signifi- 
cance may not be attached to the difference of phases unitil further confirmed. 
Moreover the value of EF, in the E and F regions may not be strictly the same, 
and the difference may be 10 or 20 per cent. To find the change of characteristics 
of daily variations with the anticipated phase-shift of EZ, and hence of W, Cases 
10 and 11 (shown in Table 4) were calculated for the phases of W (125.37) and 
(134 .37). 

As before, in these cases two maxima of v,, appear in the fore- and after- 
noon. The forenoon maximum is of an order of magnitude almost comparable 
with the afternoon maximum for Case 10 and the former is greater than the latter 
for Case 11. Moreover for Case 11 after sunset rapid increase of h,, appears for 
about one hour and is followed by the gradual decrease. These characteristics 
in the daily variations of v,, and h,, are very similar to those observed at sp-max. 
at Eq. stations. 

The course of the electron density contours vs. altitude throughout the 24 
hours period for Cases 4 and 10 are illustrated in Fig. 11. 
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(iv) To examine the conditions leading to the sunset increase of h,, which 
is a remarkable characteristic of equatorial h,, variation, Case 12 was calculated 
under the assumption of further phase lag of W. In this case the sunset increase 
of h,, is further magnified, but the time of maximum h,, lags by about 2 hours 
beyond that observed, and in the day time the variation of h,, does not resemble 
those observed at Eq. stations. One maximum of VY», appears at about 11.5hr. 

(vii) For Cases 1 to 4 the phase of W is (11>.3}) and the daily variation of 
the calculated F2 region is very similar to those observed under summer condi- 
tions in middie latitudes, and for Case 6 the phase of W is (185.0) and the daily 
variations are very similar to those observed under winter conditions in middle 
latitudes, if semi-diurnal variations are not taken into consideration. These 
characteristics will hold for the (b) type of ion production distribution as well. 
The seasonal variation of the F2 region in middle and high latitudes might thus 
be produced, if the diurnal drift velocity changes with season and the ion produc- 
tion takes appropriate forms. 

6. Discussion of the Results 

(i) When Weiss [3] discussed the effect of the semi-diurnal drift upon the 
F2 region, the amplitude of the drift was taken to be 0.59(F2)/hr., and the maxi- 
mum departure of v,, from the Chapman norm was 0.14. Here we used the drift 
velocity, which is almost an diurnal term and the amplitude of which is 0.893 $(F 2) 
—1.785(F2)/hr. Hence departures from the static Chapman norm are much 
greater than those obtained by Weiss. 

(ii) To discuss the variation of h,, Martyn [1] neglects the first and second 
terms of the right hand side of (12) and uses the equation 0N/dt=—0(N: W)/Oh. 
Therefore if use is made of W=W,sin(vi+0,,) which is independent of altitude, the 
following equation is obtained: 

Ahy, =(W,,/nw) sin (nd+9,—7/2), (18) 
wnere 0=7.3x10™°. 

If, as in the present treatment, W,=1.2410%cm/sec. is taken, the ampltude 
of 4h,, equals W,/o~170km. On the other hand, according to calculated results, 
the amplitude of 4h,, is about 60km, and is much less than that calculated by (18). 
Thus, if one uses (18) to obtain the drift velocity from the observed data, fairly 
great error will be introduced. 

(iii) According to observed results near the equator, the increase of electron 
density at dawn begins at about 5hr., but the calculated increase lags about one 
hour or more in any case, than those observed. This is due to the fact that ion 
production before ground sunrise is very small, even for type (b) in section 4. 
The early increase may be due to photo detachment by visible solar radiation 
from the negative ions formed gradually at night by molecular oxygen which is 
suggesed by Nicolet and Mange [30] to have a fairly great concentration in the F2 


region, 
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(iv) The distribution with altitude of recombination coefficient used in section 
4 may be somewhat different from the real one, so that the present calculation 
shows only a rough order of magnitude. 

(v) The electrical conductivity of the E region and its variation mentioned 
in section 3 seem to be plausible to fit the estimation by radio observation [31] 
but observations of the electron density of the E region at night are still too 
meagre to confirm this situation. 

For the electric field E, obtained by H. Maeda [20], the diurnal term is about 
seven times that of the semi-diurnal one. This result is quite different from 
Martyn’s conclusion [1] that the main term of the electric field for S, is semi- 
diurnal. The variations of v,, and h,, obtained here based on (8) agree with 
observed results. On the other hand the results calculated by Weiss [3] based on 
the semi-diurnal drift do not resemble those observed on the equator. The effect 
of the daily variation of the temperature upon the F2 region has been examined 
by Weiss [32] and Shimazaki [33] and it is seen that the effect is very small com- 
pared with those obtained here. 

7. Distribution of the #2 Region near the Magnetic Equator 

In the following, we shall evaluate the effect of the drift in the neighbour- 
hood of the magnetic euator. The main magnetic field has been taken as that of 
a dipole with its axis coincident with that of the earth’s rotation. The vertical 
drift velocity W in the F region can be written, using Equation (2) in II], as shown 
by Martyn [1], 

W =(E,/H)cos ¢= {(¥Y—uH.)/H}cos¢ 
=(Y/H)cos $+(1/2)+a-sin2¢, (19) 


where Y is the eastward component of the electrostatic field, and wz is the south- 
ward component of horizontal air velocity and H, is vertical component of H 
measured positive upwards. Below an altitude of about 140km the value W/E, 
decreases rapidly downwards except close to the magnetic equator. 

From the theory of potential, the diurnal electrostatic field expressed by 
wsin(A+96,) in (7) may be fairly uniform between +20° of latitude. However, 
Ampere force will accelerate the southward component of the air motion in the 
F2 region in such a way that the value of the diurnal component of E, and hence 
of W at about 20° is much smaller than that at the equator and the phase of the 
former proceeds that of the latter by an amount not far from 90° in a similar 
way as shown in III (§ 10). 

Referring to Case 5 of Fig. 12, a decrease of v, in the F2 region from »%,, 
at noon by the drift at 20° will be very small, and moreover the ion production in 
these latitudes is roughly the same as at the equator. According to Weiss, the effect 
of the semi-diurnal drift upon »,, is very small. It would therefore be possible 
that N,,(F2) at 20° would be greater than that at the equator. 

Complete treatment of this problem is, however, very difficult until the dis- 


Geomagnetic Distortion of the F2 Region on the Magnetic Equator 143 


tribution of electric field and air motion in the F2 region due to various causes 
is known for the whole earth. 
8. Concluding Remarks 

From the foregoing analysis the following conclusion may be drawn: 

The eastward component of the electric field producing the S, current on 
the magnetic equator produces the vertical drift of the F2 region at the same 
time, the main term of the drift velocity is diurnal. The calculated daily varia- 
tions of the F2 region, under consideration of this drift, have a striking resem- 
blance with those observed near the magnetic equator. Moreover when the 
ion production take its maximum value at about 200km, there appears a lower 
secondary maximum of electron density which agrees well with the observed F1 
layer in its density and altitude characteristics. 

The change of characteristics of the daily variation of the F2 region with 
an sp-cycle seems to be accounted for by a variation in ion production and by a 
slight shift in the phase of the diurnal drift. 

The daily F2 variations of the equatorial type have a resemblance with 
those under summer conditions in middle and high latitudes and when the drift 
velocity has such a phase (18.07) as used in Case 6, daily variations are very 
similar to the observed diurnal variation in winter conditions in the latitudes. 
Thus the seasonal variation of the F2 region in the middle and high latitudes may 
be greatly influenced by the seasonal variation of the diurnal drift. 
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Reverse and Normal Magnetism of the Basaltic Lavas at 
Kawajiri-misaki, Japan 


By Eizo ASAMI 


Physical Institute, Faculty of Arts and Science, Yamaguchi University 
(Read Oct. 31, 1953 and Nov. 3, 1954) 


Abstract 


The early Pleistocene basaltic lavas at Kawajiri-misaki (Cape 
Kawajiri) have in general reverse natural remanent magnetism 
(N.R.M.), while at a part of the place both specimens of reverse 
and normal N.R.M. are found with their positions as intermixed, 
despite that they have been taken from the same rock block, and 
with their intensities ranging from 10~‘ to 10-°c.g.s.e.m.u./g. This fact 
seems to suggest that there is much room for the discussion regard- 
ing the idea that the reverse N.R.M. of the lavas would have been 
caused by a reverse geomagnetic field assumed to have occurred at 
the age of the ejection of the lavas. 

The reverse N.R.M. of the basaltic lavas at Kawajiri-misaki situated at the 
Japan Sea coast of Yamaguchi Prefecture in west Japan was first reported by the 
present author [1] in 1954. From the thermo-magnetic studies of reversely magne- 
tized igneous rocks in Japan (including the Kawajiri-misaki basalts), T. Nagata and 
others [2] have recently inferred that the direction of geomagnetic field at the time 
of formation of some of these rocks was reverse to that of the present. 

However, the present author has left room for doubt of the above inference 
and for the purpose of further clarification of the fact he has again made a trip 
to Kawajiri-misaki and sampled more numerous basalt specimens from there. And 
the N.R.M. of these specimens have been measured in his laboratory. The results 
will be shown below. 

The age of eruption of the basaltic lavas of Kawajiri-misaki is generally 
considered to be the early Pleistocene and they are considered also generally to 
have undergone no conspicuous dynamical disturbances since their eruption. In 
Fig. 1, the marks of cross (x) represent the sampling spots where crop out solid 
blocks of lava flow and the numerals nearby represent the specimen numbers. 
The specimens were sampled from various spots as much as possible. However, 
the inner part of the cape is covered with ordinary soils and the specimens could 
not be sampled from there. Thus, the author has measured the directions and 
intensities of magnetization of these 193 specimens collected solely from the areas 
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Fig. 1 Sampling spots at Kawajiri-misaki. 
Thick numerals represent the specimen numbers of normal N.R.M. and thinn those 
of reverse N.R.M. 


fringing the coast by means of an astatic magnetometer and obtained the following 
results. 

In the first place, it is worth noting that the directions of the N.R.M of the 
specimens show a remarkable distribution; that is, in the belt C (Fig. 1) occupying 
the eastern coast and also in the belt A occupying the southern half part of the 
western coast the magnetizations are all reverse, while in the belt B occupying the 
northern half part of the western coast the normal and the reverse magnetizations 


are intermixed in positions and in the partial belt b occupying the middle part of 
the belt B the magnetizations are all normal. 
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N) 
Fig. 2 Directions of N.R.M. of the specimens in the belts A and C. 


In Fig. 2 are shown the Schmidt’s projections of the reverse magnetizations 
in the belts C and A, and in Fig. 3 those of the normal and the reverse magneti- 
zations in the belt B with the partial belt b inclusive. In the both figures the 
vector of the present geomagnetic field whose dip is +49° downward is represented 
by a thick cross (x). : 

The intensities of magnetization of the specimens are shown in Figs. 4, 5 and 
6 by the lengths of the thick horizontal lines which are arranged on the ordinate 
in order of the positions of the spots the specimens have been taken from, the 
numerals nearby representing the specimen numbers. 

The arithmetical average values of the directions and intensities of the N.R.M. 
and their ranges are shown in Table I. For the sake of reference, the similar 
quantities of the specimens from two places, one Tsuo (Fig. 1) 6km. to the east 


and the other Tawarashima 5km. to the south, of Kawajiri-misaki, are shown in 


Table II. 
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Lower hemisphere 


Upper hemisphere 


S 
Fig. 3. Directions of N.R.M. of the specimens in the belt B with the partial belt b inclusive. 


Table I N.R.M. of the Kawajiri-misaki basalts. 


Normal | Total number | Declination | Dip | (10-3 ie /g.) 
Locality or of — --—— "ae 
: Average \Average Average 
Reverse the specimens | value | Range “Value Range vatar | Range 
Belt B Normal 38 N 21°E |N 54°W +46° |4+ 9°~+80°) 5.5 0.2~28.1 
~N 67°E 
nf Reverse 36 N169°E |N124°E —57° |—38°~—75°} 1.1 0.3~ 6.7 
~N149°W 
Belt A Reverse 44 N160°E |N126°E —47° |—23°~—71"| 1.0 0.2~ 5.0 
~N137°W 
Belt C Reverse ae N168°E |N120°E —55° |—32°~—74"| 1.2 0.2~ 9.1 
~N144°W 
Total region| Reverse 155 N166°E |N120°E —53° |—23°~—75°"| 1.1 0.2~ 9.1 
~N137°W 
Table II N.R.M. of the other basalts. 
| Tsuo Normal 6 | Naver |N 8E +35° |+28°~+52°| 1.3 | 0.5~2.1 | 
| | ans ~N 36°E é heeeun 


“Tawara- ee “ar gege WRIOEE, cL ae CEO ner Feri 
aweshima | Normal | 5 | NPE IN ep] +55° |\+49°~ +58] 14 | 0.6~2.0 | 


149 


Reverse and Normal Magnetism of the Basaltic Lavas at Kawajiri-misaki, Japan 


<--—— 


Reverse N.R.M. 


—3 ¢.g.8.¢.m.u./g, 


10 


10 


me naowvr oon aoa 


Nn 


NNNANA AN C 


gts 10 ©.g.8.¢.10.u./g. 


1 


2 


Reverse N.R.M. 
<— 
3 


SP AGSAAD oo 
aoe 
aaa 


OMDMMT™ A RH ORTH NH DCH MO AMM DAH CHN | ein 
DAMA AAA CSOT OC HBF BABE AN A 
He Be et et te! 


oan tien ianianl 


CSOHAMTOR OCH MOP 
HAR RR ANNA 
Se Set et 


—3 ¢.g.2.¢.m.u./g- 
0x10 


2 


— 


Reverse N.R.M. 


Fig. 4 Intensities of N.R.M. of 


the specimens in the belt 


A. 


AMSNAAMYN MN OD RADROHAMNSOATAMHATNORNODMYAGNRSCHK HASNT OnR SH ARHHE DRS TMA CA SHE OD 
Oo oF ae © SUS Sicr'lon a9 on e908 co MH OD COS 19 SD BS OS ae SS tS Ce 8 OS A A SAD 


—3 c.g.n.c.m u/s, 


OH x 10. 


10 


< 


Reverse N.R.M. 


Intensities of N.R.M. of the specimens 


in the belt C. 


5 


Fig. 


150 E. Asami 


Reverse N.R.M. <— — » Normal N.R.M. 


—3 c.g.e.c.m.u./g. 


10) 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 10 20 30x10 
87 
85 
20 \ 
19 
18 
: 17 
82 
1 
; 80 
79 
78 
77 
76 
75 
1 
; 74 
73 
7 
4 15 
70 
6 
F 14 
13 
68 
12 
67 
66 
65 
64 
62 
= 1l 
10 
‘61 
60 
59 
58 
57 
s 55 
54 
53 
52 
51 
50 
49 
9 
8 
7 
aes 6 
fee oy 48 
Cs {7 
46 
5 
5 
44 
4 
41 
40 
39 
38 
36 
35 
a 
2 
1 
34 
33 
32 
29 
26 
24 
23 
22 
21 


LU Ds eB Tg ee Wire: Cas ee ey | 2 8 4 8 ely Salo mpines 30 x 1973 Cee 
Reverse N.R.M. <— —- Normal N.R.M. 
Fig. 6 Intensities of N.R.M. of the Specimens in the belt B with the partial belt b inclusive. 


Reverse and Normal Magnetism of the Basaltic Lavas at Kawajiri-misaki, Japan 151 


Thus, the characteristics of the directions and intensities of the N.R.M. of 
the basaltic lavas at Kawajiri-misaki have been definitely clarified. Above all, the 
author has now recognized the decisive fact that the lavas at Kawajiri-misaki 
have in general reverse N.R.M., but at a part of the place both specimens of 
reverse and normal N.R.M. are found with their positions as intermixed, despite 
that they have been taken from the same rock block where any. geological distur- 
bances can hardly be noticed; for example, two specimens of normal No. 80 and 
reverse No. 81 have been taken from the spot considered to be the entirely same 
_ Tock block in the northern part of the belt B (the distance separating the positions 
of these two specimens is only a few meters). 

Moreover, as shown in Figs. 4, 5 and 6 and Tables I and II, the intensity of 
magnetization of these basalt specimens at Kawajiri-misaki is not uniform, while 
at all the other places [3] of Yamaguchi Prefecture, such as Tsuo, Tawarashima, 
it is fairly uniform and the N.R.M. is wholly normal. Roughly speaking, the average 
value of the intensities of magnetization of the specimens from Kawajiri-misaki is 
10-*c.g.s.e.m.u./g. for the normally magnetized basalts of the partial belt b, 10-* for 
the reversely magnetized basalts of the total region and 107? for the normally 
magnetized basalts at the spots sandwiched between the spots of the reversely 
magnetized basalts of the belt B. 

On the other hand, from the results of the thermo-magnetic study of our 
Kawajiri-misaki basalts, it has already been reported [2] that of both the reverse 
and normal specimens the N.R.M. is mainly due to the thermo-remanent magnetism 
and not to the local anomalous magnetization owing, for example, to thunderbolts. 

It will be worth while to note that these facts may have a significant mean- 
ing for the study of palaeomagnetism. Generally speaking, if the igneous rocks 
have kept their N.R.M. stable since they acquired it as the thermo-remanent 
magnetism under the effect of the geomagnetic field in the remote epoch, it might 
be inferred that the direction of the N.R.M. represent that of the geomagnetic 
field at the time of formation of the rocks. It is generally believed that the 
thermo-remanent magnetism of stable ferromagnetic minerals only is able to remain 
stable during a long time. The above-mentioned facts seem to suggest that there 
is much room for the discussion regarding the inference [2] that the reverse N.R.M. 
of the lavas at Kawajiri-misaki would have been caused by a reverse geomagnetic 
field assumed to have occurred at the age of the ejection of the lavas. Because, 
should this inference be taken for granted, how could we explain the fact that at 
a part of the place the normal and the reverse magnetizations are intermixed in 
positions in the same rock block? If the normal or reverse of the directions of 
N.R.M. of rocks depend on that of the geomagnetic field, all specimens from the 
same rock block must show without exception only one side of either normal or 
reverse N.R.M.. 

However, it would be a little questionable whether the position intermixing 
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of the reverse and the normal magnetizations has taken place in the really same 
rock block or not. A key may, therefore, be afforded by a close geological field 
observation for confirming whether the normally magnetized rocks are dyke rocks 
penetrating the body of the reversely magnetized rocks or not. A more proper 
key would be the exact age determinations of both the reversely and the normally 
magnetized rocks. 

On the other hand, another interpretation seems to be given, on the assump- 
tion that the above-mentioned fact has taken place in the really same rock block, 
by the idea most recently proposed by N. Kawai and others [4] that the magnetism 
of rocks cannot always be kept stable and the self-reversal of magnetization of 
rocks is possible to occur by a solid-phase transformation in ferromagnetic minerals 
responsible for the remanent magnetism of the rocks through a long period in 
geological time scale. The author has obtained the dependency of the intensity 
of saturation magnetization upon temperature for the several specimens from 
Kawajiri-misaki by means of a magnetic balance. From the result, it seems to be 
not particularly difficult to suppose that the facts found at Kawajiri-misaki may 
be in connection with this idea to some extent. 

In concluding, the author wishes to express his sincere thanks to Dr. M. 
Matuyama, M.J.A., President of Yamaguchi University, for his kind direction 
throughout this research. His cordial acknowledgement is also due to Prof. N. 
Kumagai and Prof. M. Hasegawa of Kyoto University and Prof. T. Nagata of Tokyo 
University for their suggestive guidances and encouragements. He is also indebted 
to Mr. H. Domen for his assistance in this work. 
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LETTERS TO THE EDITOR 
On the Solar Lyman Beta Radiation and the Ionosphere 


It has been known that there exists an atomospheric window at the wave- 
length of the Lyman beta line of hydrogen in the ionosphere [1]. R. v. d. Woolley 
suggested that the solar Lyman beta radiation emitted from the chromosphere 
would be an effective radiation which gives rise to the radio fade-out [2]. 

The Lyman beta line of hydrogen is, however, coincident with the radiation 
caused by transition between the ground 2°P level and the excited 3°D level of 
atomic oxygen in wavelength. {(Z,4 1025A, OI 2°P-3°D 41026A) As atomic oxygen 
is one of the main constituents of the ionosphere in the E and F regions, the 
incident solar Lyman beta line will be absorbed by atomic oxygen according to 

O(2°P)+ hy>O(3'D) (1) 

The produced O(3°D) atom further passes partly to the continuous level by 
absorbing the visible radiation (47813A) and partly to the 3°P level, emitting the 
infrared radiation (4 11287A) i.e. 

O(3°D)+hv(<47813A)>O*+e, (2) 
O(3°D)->0(33 P)+hv(A 11287A). (3) 
Further we have the following 
processes 
O(3°P)—>O(3'S)+hv(A 886A), (4) 
O(3*S)—>O(2°P)+hv(A1302A), (5) 
O(3°S)+hv( <2 3033A)>O*+e, (6) 
O(2°P)+hv(>13.58eV)>O*+e. (7) 

These processes are to de- 
termine the feature of the inci- 
dent solar Lyman beta radiation 
through the ionosphere. 

All of these processes ex- 
cept (2) (6) and (7) make the inci- 
dent Lyman f-line split into the 
two scattered radiations i.e. the 
radiatons 2°P-3'D and 2° P-3°S, 
accompanying absorption and 
emission of the visible as well as 
the infrared radiations. 

Here we must pay an atten- 
tion to these radiations. For the 
thermal emission at these wave- 
lengths should be taken into ac- 
Energy Levels of OI count in the F2 layer where tem- 
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perature has been estimated to be 10° °K in the order of magnitude. For instance, 
with temperature of 10°°K the ratio of the thermal emission from the ionsphere 
to the incident solar radiation par c.m’. and sec. is 0.3 at 8446A and 1.0 at 11287A; 
with 2x10°°K it becomes 5.7102 and 1.510? respectively, provided that the solar 
radiation at these wavelengths is that of black body with temperature of 6x 10? 
Therefore, in such a layer with high temperature the splitting of the incident 
radiation into two seems to be far less than in the lower layers with lower temper- 
ature. 

The processes (2) and (6) produce an ion of atomic oxygen through visible 
radiation. The produced ion will recombine with electron through the radiative 
recombination or some other recomibinaion reaction suggested by many authors. 
S. Miyamoto put forward the theory that in the ionosphere recombination radiation 
would be available for further ionization [3]. This theory has been developed later 
by Y.Inoue and proved to be relevant [4]. He has considered all possible recom- 
bination reactions from theoretical as well as observational standpoints and con- 
cluded that the most probable processes are the radiative recombination and the 
mutual-neutralization, both of which have the same effect on the radiation field 
apart from their rates.and contribute to the scattering of the ionizing as well as 
the split radiations [4]. 

We assume that the ionosphere principally consists of atomic oxygen and 
molecular nitrogen in the E and the F regions. The ionization potential of molecular 
nitrogen is 15.58 eV, while that of atomic oxygen being 13.58 eV. Furthermore 
there is no considerable absorption at wavelengths of the two split radiations. 
Molecular nitrogen, therefore, does not play any important role in our problem. As 
being the case, the radiation field has no sink in hydrodynamic term throughout 
the ionosphere down to the region where the appreciable absorption of molecular 
oxygen takes place. 

Now we examine the splitting of the incoming solar Zf radiation in the lower 
region following the theory of radiative transfer. 

It is assumed that the radiation field is in a quasi-stationary state and varies 
with the zenith distace of the sun in a slow manner. The equations of radiative 


transfer are as follows, 
cos0 He — wa( Te NA ), 
dt 


N,Bip 2 
dI,, NsA 
0 ere : geen PS44S1 
cos dr te (1. N.Bis ), 
dl, NcN,A 
p22o — yp, 4VCtVasicr 
$98 dt Ic N,Byic 


where Bip, Bis, Bic, An1, Asi and Ac; are Einstein coefficents of neutral oxygen 
between the two levels indicated by the suffixes. We denote the four levels i.e. 
2°P 3°D 3°S and the contionus levels by 1, D, S, and C, respectively. NV, is electron 
density, and N,, N,, Ns and Nc are number densites in the respective levels. The 
radiation between 1 and D is named a- radiation, between 1 and S a'-, and between 
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1 and D C-radiation. + is optical depth and is defined by 
dr= —N,B,cedz, 
where z is a coordinate taken vertically upwards. 


w, = Pines 
I Bice 
ens Byis€ a: 
ap emit: A 
ic€é 


where €,, €,, and € photon energy of a-, a’- and C-radiations respectively. As to 0 and 
I, we follow the usual definition in the theory of radiative transfer in case of the 
plane parallel layer atomosphere. 

The population of each level is determined by the following cyclic equations 


me =NiBirJat+NeNeAcot+NsBs x] 9 —No(Anit Boclp + As); 
oer =N,BysJa+NceN.Acst+tNpAps —Ns(Asit Bsolor+Bs oy), 
oa =NpAptNsAsitNcN Aci—Ni( Bidet Bisa t+ Bile); 
oe =N,Bicle+NoBoc)s+NsBscJp:—NeN(Acit Acot Aes), 


where the subscripts §, 8’ and f” pertain to the visible radations between D-C, 
S-C and S-D respectively. We take four levels, 1, S, D, and C levels into con- 
sideration although there is one more level 3°P located between S and Dlevels. This 
procedure would be valid approximately, because thermal radiation from the iono- 
sphere plays a role only in the F2 layer and its eftect rapidly vanishes in the lower 
region. And all visible as well as infrared radiations concerning our problem come 
from outside. We assume their effects to be constant on account of their abundant 
quantities. Then, the 3*P level does not appear explicitly, but implicitly through 
Aps and Bs, in the equation. We define as in usual case J., Jar Jc, Js, Jgr and Jan 


by 
1 ees! 
Ja= ig |e Je= a Id, etc. 


J. is composed of the two parts, the incident part J%.exp(—rtw.sec¢) and the scatter- 


ed part J, as follows 
Je=Jéaexp(—tw.sect)+/, 
where J%. is the solar radiation at r=0 and ¢ is the zenith distance of the sun. 
The Edington approximation is asumed to be applicable to the scatted part. 
Finally we get for the region with r>1 which corresponds to the EF and Fl 


layers 


AatAcn _ , Bsc: re BgpJ an 
Ja _ War - AcitAcnotAcs Asi Asi 
Jat. Wy Amt+Acs BocJs br Aps 


AcgitAcnt+Acs 1 Ass An 


Although the coefficients in this formula have not yet been known exactly, the 


following two cases would be possible. 
Case 1 The f#” radiation is more effective than the @ and ’ radiation. 
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Jo ~ Bas é Ap ; BspJ yx lO ee 
Ja ——Ass By Ans 


Case 2 The processes with # and /’ radiations play an important role. 
w- WA gee eee o 
pa 

From these results we conclude that the incident solar Lyman beta radiaton 
is almost changed into a’-radiation in this region. 

Recently the first successful observation of the Lyman alpha radiation was 
reported by W.B. Pietenpol and others [5], [6]. Once such an observation is made 
at higher altitude and in shorter wavelength region our estimation will directly be 
checked. It is interesting to suggest in this connection that the Lyman beta radi- 
ation would hardly be observed except in the F2 region with high temperature 
and further whether or not the a’-radiation with a little longer wavelength than 
Lyman alpha is observed should be a critical test for our scattering model of the 


ionosphere. 


The author wishes to express his sincere thanks to Profs. M. Hasegawa and 
S. Miyamoto for their kind direction and encouragement through the work. The 
author is also indebted to Mr. Y. Inoue for his invaluable suggestions and discus- 
sions. 
Susumu Kato 


Geophysical Institute, Kyoto University, 
Kyoto, Japan 
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On Latitudinal Distributions of Diurnal and Semi-diurnal 
Components of h’F2 of the Ionosphere 


A study has been made of the latitudinal distributions of the phase and 
amplitude of the diurnal and semi-diurnal components of h/F2, the minimum virtual 
height of the F2 region. 

The data subjected to analysis are the hourly median values of h’F2 from 36 
world-wide ionospheric observatories in December 1945, and from 39 in June 1946. 
A harmonic analysis of #’F2 in Equinox has already been made by Shimazaki [1], 
so that we ommitted the analysis of the data in this season. 

The results of the a analysis are shown in Figs. 1-8, in which the components 


are distributed against the geographic and geomagnetic latitudes. General charact- 
eristics are summarized as follows: 


(i) Diurnal Component 
(a) The amplitude is larger in northern hemisphere than that in southern 


in summer and vice versa in winter. However, the very large value in high latitude is 
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Fig. 1 Distribution of amplitude of h'F2 Fig. 2 Distribution of amplitude of h'F2 
diurnal component in December diurnal component in June 1946, 
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Fig. 3 Distribution of phase of h’F2 diurnal 
component in December 1945, plotted 
against geogarphic and geomagnetic 
latitudes. 
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Fig. 5 Distribution of amplitude of h’F2 
semi-diurnal component in December 


1945, plotted against geographic and 
geomagnetic latitudes. 
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Fig. 4 Distribution of phase of A’F2 diurnal 
component in June 1946, poltted against 
geographic and geomagnetic latitudes. 
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considered to be unreliable, judging from the behaviour of the height of the maxi- 
mum electron density. 

(b) The time when the positive maximum occurs is nearly noon in northern 
hemisphere and midnight in the southern in summer and vice versa in winter. In 
geomagnetic equator, it seems that there is time lag of 4-6 hours. 

(ii) Semi-diurnal Component 

(a) The distribution of the amplitude of this component through two hemis- 
pheres is analogous to that of the diurnal component, but the magnitude of the 
former is fairly small compared with that of the latter. 

(b) The time when the positive maximum occurs is nearly midnight and noon 
in summer and winter through two hemispheres. It is considered that, adding the 
result of analysis by shimazaki, the phase is constant through the year world- 
widely. There are in some stations at high latitude where the phase lags about 


4-6 hours. 
By Teruo SATO and Tomikazu NAMIKAWA 


Geophysical Institute, Kyoto University 
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